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1 .  Introduction 

We  consider  the  problem  of  allocating  a  number  of  identical  resources  to  requests  arriving  at  the 
sites  of  a  distributed  network.  We  assume  that  the  network  is  configured  as  a  tree.  The  ncdes  of  the 
tree  are  processors  and  the  edges  are  communication  lines  connecting  the  processors.  Processes  at 
a  node  may  communicate  only  over  the  tree  edges,  with  processes  at  other  nodes.  Resource 
allocation  is  managed  by  a  collection  of  communicating  resource  allocation  processes,  one  at  each 
node.  We  will  henceforth  refer  only  to  the  node,  identifying  it  with  both  the  processor  and  the 
resource  allocation  process  at  the  node. 

From  time  to  time,  a  request  arrives  at  a  node  (potentially  any  node  of  the  network)  from  the 
outside  world.  One  of  the  resources  should  eventually  be  granted  to  the  request,  subject  to  the 
following  conditions: 

1.  No  resource  is  granted  more  than  once.  (Once  granted,  a  resource  is  not  returned.  Thus,  the  ^ 
is  no  legitimate  reason  to  grant  it  more  than  once.) 

2.  At  most  one  resource  is  granted  to  each  request. 

3.  A  node  grants  resources  only  to  those  requests  which  arrive  at  that  node. 

4.  If  the  number  of  requests  is  no  greater  than  the  number  of  resources,  then  each  request 
eventually  receives  a  resource. 

5.  If  the  number  of  resources  is  no  greater  than  the  number  of  requests,  then  each  resource  is 
eventually  granted  to  a  request. 

For  convenience  in  describing  allocation  of  specific  resources  to  specific  requests,  we  assume 
that  each  resource  and  request  has  a  unique  identifier. 

The  execution  model  for  this  distributed  network  is  event-based.  Two  types  of  events  may  occur 
at  a  node:  (1)  a  request  may  arrive  Irom  the  outside  world,  and  (2)  a  message  may  arrive  from  a 
neighbor  in  the  tree.  Each  event  triggers  an  indivisible  step  at  the  node.  This  step  may  include 
changing  state,  sending  messages  to  other  nodes,  and  granting  resources  to  requests.  (We  ignore 
the  time  involved  in  this  local  processing  when  we  measure  the  response  time,  considering  only  the 
communication  time.)  We  assume  that  the  communication  lines  are  reliable,  that  is,  each  message  is 
delivered  exactly  once.  Hov/ever,  we  do  not  make  any  assumptions  about  the  order  ot  message 
arrivals. 

There  are  many  interesting  approaches  to  solving  this  resource  allocation  problem.  In  a 
centralized  approach,  all  resources  are  controlled  by  a  single  central  node.  When  a  requosi  arrives  at 
a  possibly  different  node,  a  "buyer"  is  commi.ssioned,  who  travels,  via  messages,  to  the  ce  itral  node 
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to  obtain  a  resource.  The  buyer  then  carries  the  resource  back  to  the  node  where  the  request 
originated,  so  that  the  resource  can  be  granted  to  the  request  at  that  location. 

An  alternative  approach  is  to  decentralize  control  of  the  resources,  giving  each  node  of  the 
network  control  of  some  of  them.  In  this  approach,  the  buyers  must  search  for  the  resources.  An 
important  choice  to  be  made  in  designing  an  efficient  search  strategy  is  the  choice  between  sending 
only  one  buyer  to  search  for  resources  for  each  request  and  sending  several  buyers  in  parallel  to 
search  different  parts  of  the  tree.  The  former  search  strategy,  which  we  call  the  sequential  search 
strategy,  avoids  a  number  of  problems  arising  from  the  parallel  strategy,  such  as  what  to  do  about 
other  buyers  when  one  of  them  has  found  a  resource.  The  next  choice,  if  the  sequential  search 
strategy  is  used,  is  the  choice  of  direction  to  search  the  tree.  A  good  choice  would  involve  guessing 
which  nodes  are  most  likely  to  have  free  resources  when  the  buyer  arrives  at  them. 

ether  strategies  involve  combining  a  decentralized  search  strategy  with  a  dynamic  resource 
redistribution  str  it  xjy,  letting  resources  search  for  requests  (rather  than  vice  versa),  or  giving  nodes 
control  o!  fractions  of  rucources  rather  than  whole  resources. 

One  cornpl'''/ity  me. '.sure  wh'ch  is  useful  for  evaluating  different  strategies  is  the  expected 
response  time.  This  is  a  measure  upon  vuiich  any  of  the  design  choices  could  have  a  major  impact. 
For  example,  the  response  time  when  using  a  centralized  strategy  must  depend  strongly  on  the 
network  size  Hov/ever  the  decentralized  strategies  have  the  potential  of  depending  on  this  size  to  a 
lesser  extent. 

In  the  first  half  of  tnis  paper,  we  present  an  algorithm  for  solving  this  resource  allocation  problem. 
Our  ..ilqorithm  is  a  decentralized  solution  in  which  each  node  controls  some  whole  number  of 
resources,  A  sequential  search  strategy  is  used,  in  which  the  direction  to  be  searched  is  chosen 
riondeterminiulically.  Certain  neighborhoods  of  the  node  at  which  a  request  originates  are  exhausted 
before  the  s^-arch  proceeds  to  mote  distant  neighborhoods. 

In  erdef  to  gain  s.ome  insight  into  the  expected  respoaso  time  for  our  algorithm,  v.o  simulated  its 
behavior  in  some  sp’.cial  cases  Ttia  nondeterministic  choice  of  search  direction  was  resolved  by 
u,sin'j  a  probaiiiiisac  cnon,','.  '.vhere  the  probabilities  for  the  different  directions  depended  on  the  initial 
olaceuient  of  resourc'.s  in  these  directions.  We  assumed  an  exponential  distribution  for  time  of 
atr,v;J  of  rofjur'sts.  a  uniform  distribution  for  arrival  location,  and  a  normal  probability  distribution  for 
m  -S'-  1  i  I  .'(..iy  time,  'vVe  also  assumed  that  all  leaves  of  the  network  tree  were  at  tlie  same  distance 

i.  on;  tliu  root.  amJ  that  ro'Cjueids  and  msources  occurred  only  at  le.'ives.  We  first  noted  that  expected 
r•,‘^.pon5•;  lime  was  exiremely  good,  with  an  upper  bound  that  seemed  to  be  independent  of  the  size  of 
the  no-t'.vmk,  this  w.is  m  m.arkod  contra.st  to  a  ce-ntrali/od  algorithm.  Next,  we  made  a  surprising 
ot.s‘1  V ah  ov  the  oxpr-ch  d  looponso  time  appeared  to  be  a  nondecroasing  function  of  the  expected 
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interarrival  time  for  requests.  If  true,  this  observation  would  imply  that  the  worst  case  for  the 
algorithm  was  actually  the  case  where  requests  come  in  so  far  apart  that  they  are  processed  one  at  a 
time.  This  observation  contradicted  our  preliminary  intuitions  about  the  algorithm:  we  had  thought 
that  the  worst  cases  would  arise  when  there  was  greatest  competition  among  requests  searching  for 
resources. 

Using  these  observations  as  hints,  we  were  able  to  carry  out  a  substantial  amount  of  analysis  of 
the  algorithm's  behavior,  and  this  analysis  comprises  the  second  half  of  this  paper.  Namely,  we  prove 
an  upper  bound  on  the  expected  response  time  for  a  special  case  in  which,  among  other  restrictions, 
all  leaves  o*  the  network  tree  are  at  the  same  distance  from  the  root,  and  requests  and  resources 
occur  only  at  leaves.  First,  we  show  that  the  response  time  can  be  bounded  in  terms  of  the  number  of 
messages  cf  one  type  that  are  sent  during  the  execution  of  the  algorithm.  Then  we  show  that  this 
number  of  messages  is  a  nc  decreasing  function  of  the  interarrival  time  for  requests.  The'efore,  the 
worst  case  occurs  v.hen  requests  come  in  so  far  apart  that  they  are  processed  sequentially.  We 
analyr-.-  the  expected  time  for  the  sequential  case,  showing  it  to  be  bounded  by  a  constant, 
independent  of  the  s'ze  of  the  netv;ork.  It  follows  that  the  expected  response  time  for  the  a  gorithm  is 
also  boLindeo  by  a  constant. 

Although  the  e,»pected  response  time  for  our  algorithm  is  very  good,  we  do  not  claim  that  it  is 
optimal.  In  fact,  there  are  some  simple  changes  that  one  would  expect  to  yield  improvements. 
Unfortunately,  with  these  changes,  the  algorithm  can  no  longer  be  analyzed  using  the  same 
techniques;  thus,  v.'o  are  not  really  certain  that  they  are  improvements  at  all. 

There  are  several  contributions  in  this  paper.  First,  we  think  that  the  algorithm  itself  is  interesting. 
Second,  we  have  identifier!  an  interesting  criterion  for  the  performance  of  a  distributed  algorithm: 
th.it  the  pi-rformance  be  iriclct.'cndcnt  of  the  size  of  the  network.  Satisfying  this  criterion  seems  to 
retjuiro  an  appropriate  decentralizerl  style  of  programming.  Third,  the  analysis  is  decomposed  in  an 
inh:ri’stirKj  v/ay:  a  '-eciuonli.il  version  is  analyzed  using  traditional  methods,  and  the  performance  of 
the  concurrent  nlijoi  ithm  is  sho'.vn  to  be  bounded  in  terms  of  the  sequential  algorithm.  It  is  likely  that 
this  kind  of  decomposition  will  prove  to  be  useful  for  analysis  of  other  distributed  algori.hms.  For 
instanr;e.  a  similar  ileconiposilion  w;is  used  in  the  proof  of  correctness  of  a  systolic  stack  [Guibas, 
an(  I  I  ring  ( !  9b?)]. 

Fho  contents  of  the  rest  of  Itm  papirr  are  as  follows.  Section  2  contains  the  algorithm,  and  Section 
;3  coiitams  argument:;  for  ds  corroctn<'ss.  Sections  4-6  contain  the  analysis  of  the  algorithi  i.  Section 
4  proves  Ihe  monotenic  ity  result,  wliich  implies  that  the  sequential  case  of  tlie  algorithm  is  worst. 
•S'  ‘(.linn  5  analyze'.  |h  '  ■  rpjr'otial  r.ase.  Section  6  pulls  together  the  results  of  Sections  4  and  5,  thus 
giving  ;i  (jeneral  upp.  t  bound.  Fiis'illy,  Section  7  descrities  some  remaining  questions. 


2.  The  Algorithm 

In  this  section,  we  present  our  algorithm.  We  begin  with  an  informal  decription,  followed  by  a 
more  formal  presentation. 

2.1.  Informal  Description 

We  assume  that  the  network  is  a  rooted  tree. 

Our  algorithm  is  a  decentralized  algorithm  with  a  sequential  searching  strategy.  Requests  send 
buyers  to  search  for  resources.  When  a  buyer  finds  a  resource,  it  "captures"  it.  Each  captured 
resource  traveis  back  to  the  origin  of  this  buyer  (or  possibly  some  other  buyer,  if  there  is  interference 
between  the  processing  of  concurrent  requests),  so  that  the  grant  can  occur  where  the  request 
originated. 

When  a  request  or  buyer  arrives  at  any  node,  any  free  resource  at  the  node  is  captured.  If  there 
are  no  free  resources  there,  a  buyer  is  sent  to  a  neighboring  node,  determined  as  follows.  Each  node 
keeps  track  of  the  latest  estimate  it  knows,  for  the  number  of  resources  remaining  in  each  of  its 
subtrees  Each  node  sends  a  message  informing  its  parent  of  each  new  request  which  has  originated 
within  the  child’s  subtree.  The  estimate  'A'hich  a  node  keeps  for  the  number  of  resources  remaining  in 
a  subtree,  is  calculated  from  the  initial  placement  of  resources  in  that  subtree,  the  number  of  requests 
v;hich  are  known  to  have  originated  within  that  subtree,  and  the  number  of  buyers  which  the  node  has 
already  sent  into  that  subtree.  In  order  to  decide  on  the  direction  in  which  to  send  a  buyer,  a  node 
uses  the  follov/ing  rules.  First,  it  never  sends  a  buyer  out  of  its  subtree  if  it  estimates  that  its  subtree 
still  contains  a  resource  Second,  it  only  sends  a  buyer  downward  to  a  child  if  it  estimates  that  the 
child's  subtree  contains  a  resource.  Third,  if  there  is  a  choice  of  child  to  which  to  send  the  buyer,  the 
node  maKcs  a  nondeterministic  choice.  (Later,  we  will  constrain  this  decision  to  a  probabilistic 
choice  using  a  particular  random  choice  function.  This  constraint  will  be  important  for  the  complexity 
analysis,  but  is  not  needed  for  the  correctness  of  the  algoiilhm.) 

it  is  easy  to  see  that  any  subtree’  which  a  node  considers  to  contain  no  resources,  actually 
cont.iins  no  r  jsisurces  riius  no  buyer  is  ever  sent  out  of  a  subtree  actually  containing  a  resource. 
On  the-  othi;r  hand,  the  perceived  information  about  the  availability  of  a  resource  in  a  child's  subtree 
can  be  an  uverestnnate,  in  case  of  interference  among  concurrent  requests. 

EXAMPLE 

Suppose  that  request  A  enters  at  the  node  shown  below,  and  its  buyer  travels  upward  until  it 
reaches  an  ancestor  that  perceives  the  av.iilability  of  a  resource  in  one  of  its  subtrees.  Tften  the 
buyer  travels  rlowovyarc!  tow.ard  ttiat  resource.  Shortly  heforo  A's  buyer  reaches  the  resource, 
.inotlKii  refjii'.'  ,t  (I,  .irrives  at  the  node  .sliown.  Suppose  O's  buyer  reaches  the  n.soureo  and  <  .rptuies 


it  before  A  s  buyer  does.  When  (or  before)  A's  buyer  finally  arrives  at  the  resource's  location,  it  will 
encounter  the  information  that  the  resource  is  no  longer  there.  Then  A's  buyer  will  be  sent  upward, 
backtracking  in  its  search  for  a  resource. 


Figure  1 

Although  such  interference  can  cause  backtracking,  the  buyer  will  eventually  find  a  ressource  if 
one  exists.  This  is  because  no  buyer  ever  leaves  a  subtree  actually  containing  a  resource. 

Several  optimizations  are  incorporated  into  the  algorithm,  as  follows. 

1 .  Buyers,  unlike  requests,  need  not  be  uniquely  identified.  Instead,  each  node  keeps  track  of  the 
number  of  buyers  received  and  sent  and  the  net  flow  of  buyers  over  each  of  its  incident  edges. 
Captured  resources  then  travel  in  such  a  way  as  to  negate  net  flow  of  buyers,  and  because  a  buyer 
will  eventually  leave  a  subtree  which  does  not  contain  a  resource. 

2.  Buyers  can  travel  "discontinuously".  Assume  node  v  sends  a  buyer  to  a  child  node  w,  thinking 
that  there  is  an  available  resource  in  w’s  subtree.  Assume  that,  soon  thereafter,  v  receives  a  message 
from  w,  informing  v  of  an  arrival  of  a  new  request  in  w’s  subtree,  and  implying  that  v's  previous 
supposition  of  an  available  resource  was  false.  Then  v  knows  that  w  will  eventually  send  some  buyer 
back  up  to  V,  at  which  time  v  should  send  the  buyer  in  another  direction.  Since  v  knows  this  will 
eventually  occur,  v  need  not  actually  wait  for  the  buyer  to  arrive  from  w;  it  can  create  a  new  buyer  and 
send  it  in  anticipation  of  the  later  return  of  the  first  buyer.  Since  the  first  buyer  will  not  find  any  free 
resources  in  the  subtree,  this  extra  parallelism  does  no  harm.  In  fact,  with  this  optimization,  it  is  no 
longer  necessary  for  w  to  return  the  buyer  at  all,  since  v  must  ignore  it  when  it  returns  to  it  in  any  case. 

3.  If  each  node  knows  how  many  resources  were  initially  placed  in  each  of  its  children's  subtrees, 
then  it  is  not  even  necessary  for  explicit  buyers  to  be  sent  upward  at  alll  All  that  is  necessary  is  for 
nodes  to  send  "ARRIVAL"  messages  upward  to  their  parents,  informing  them  of  the  arrival  of  new 
requests  in  their  subtree.  The  parent  is  able  to  deduce  the  number  of  resources  which  the  child 
would  like  to  have  sent  down  (i.e.  the  number  of  buyers  emanating  from  the  child’s  subtree)  from  the 
initial  number  of  resources  in  the  subtree,  the  number  of  arrivals  in  the  subtree  and  the  number  of 
buyers  already  sent  down  into  the  subtree.  We  will  say  more  in  a  moment  about  how  this  deduction  is 
made. 
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If  information  about  newly-arrived  requests  (in  the  form  of  "ARRIVAL"  messages)  only  flows  upward 
in  the  tree,  there  is  no  way  that  a  child  can  deduce  that  its  parent  would  like  it  to  send  a  resource 
upward.  Thus,  it  is  still  necessary  to  send  explicit  buyers  downward.  Let  us  designate  these  explicit 
downward  buyer  messages  as  "BUYER"  messages.  Thus,  the  algorithm  only  uses  two  kinds  of 
messages  to  search  for  resources:  "ARRIVAL"  messages  flow  upward  to  inform  parents  about  new 
requests,  and  'BUYER"  messages  flow  downward  to  inform  a  child  that  its  parent  would  like  the  child 
to  send  up  a  resource 

The  precise  deduction  which  a  parent  can  make  about  the  number  of  buyers  emanating  from  a 
child's  subtree  is  as  follows. 

Let  a  be  the  number  of  "ARRIVAL"  messages  which  have  been  received  by  the  parent  from  the 
child  Let  b  be  the  number  of  "BUYER"  messages  which  have  been  sent  by  the  parent  to  the  child. 
Let  p  be  the  number  of  resources  initially  placed  in  the  child’s  subtree.  Then  the  number  of  buyers 
perceived  as  emanuhng  from  a  chad's  subtree  is  max(a  -r  b  -  p,  0).  This  number  is  called  the  estimate 
of  "viitual  buyers  "  e.r,anating  from  the  subtree. 

That  IS.  if  th:,'  mtal  number  of  ARRIVAL"  and  "BUYER"  messages  indicated  above  is  no  greater 
than  the  initial  placement,  no  frjyers  are  perceived  as  emanating  from  the  subtree.  On  the  other 
hand,  if  this  total  is  greater,  then  tno  excess  is  perceived  to  be  the  number  of  buyers. 

Analogously,  the  riiild  node  deduces  an  estimate  of  the  number  of  "virtual  buyers"  it  has  sent  out 
of  its  siiiitree,  as  follows.  Let  a  bo  the  number  of  "ARRIVAL"  messages  which  the  child  has  sent  to  its 
pari.mt  l.i.R  b  be  tiie  number  of  "BUYER"  messages  which  have  been  received  by  the  child  from  its 
p.ar'int  I  et  p  be  the  number  of  resources  initially  placed  in  the  child's  subtree.  Then  the  number  of 
Lu.ors  lie  child  perceives  that  it  hd.s  sent  out  of  its  subtree  (also  called  the  estimate  of  "virtual 
ii'.iyers  "  .sent  out  of  the  subtree)  =  n.a<{a  »■  b  •  p.  0).  Because  of  message  delays,  the  child  and  the 
p-'.ront  may  d'ffer  on  tlioir  estin'atos  of  the  number  of  virtual  buyers. 

in  order  tc  maxo  the  actual  grants  to  specific  requests,  it  seems  necessary  that  each  specific 
identifiiyi  resource  "travel"  to  a  pom!  of  request  origin,  in  order  to  get  properly  paired  with  a  request. 
Ttua  travel  reouircs  a  tlnrrj  kind  of  message  to  be  sent  around,  namely,  a  specific  "captured 
r.-s  ,'Urc.j".  Ffie  algorithm  which  scuds  resources  around  is  particularly  simple  ■  resources  are  just 
a  .:  t  in  '  Lich  a  v.'iiy  as  to  ncgalM  the  net  flow  of  buyers.  This  |)art  of  the  algorithm  executes 
:onr  ui  v.'iU-!  hut  h,:i ;  no  effect  on,  the  searcfiing  pait. 
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2.2.  Formal  Description 

In  this  subsection,  we  present  a  program  implementing  the  algorithm  described  above.  A  sketch 
of  a  correctness  proof  is  presented  in  the  next  section.  Primarily,  the  proof  consists  of  showing  the 
correctness  of  the  invariant  assertions  made  at  various  points  in  the  program.  The  reader  rray  wish  to 
examine  the  proof  while  reading  the  program. 

We  assume  that  the  network  is  described  by  a  rooted  tree  T.  For  uniformity,  let  the  root  of  T  have 
an  outgoing  upward  edge.  (Messages  sent  along  this  edge  will  never  be  received  by  anyone.)  We  can 
then  Vyirite  a  single  program  for  all  the  nodes  of  T,  including  the  root. 

Let  V  denote  the  set  of  vertices  of  T,  Let  RESOURCES(v)  denote  the  resources  placed  at  vertex  v, 
for  each  v  C  V,  and  let  PLACE(v)  =  |RESOURCES(v)|  for  all  v.  Let  REQUESTS(v)  denote  tha  requests 
arriving  at  v.  We  assume  that  all  the  sets  RESOURCES(v)  and  REQUESTS(v)  are  finite.  Let 
PARENT(v)  CHlLDREN(v).  DESCENDANTS(v)  and  NEIGHBORS(v)  denote  the  designated  vertices 
and  sets  of  vertices,  for  vertex  v. 

The  kinds  of  messages  used  are  "ARRIVAL",  "BUYER"  and  messages  corresponding  to  specific 
captured  resources. 

Program  for  node  v,  v  €  V: 

In  the  program  for  node  v,  we  use  RESOURCES  as  a  shorthand  for  RESOURCES(v),  and  similarly 
for  the  other  notation  above. 

It  IS  convuo  '"'It  to  think  of  the  state  of  v  as  consisting  of  "independent  variables"  and  "dependent 
variatdes".  The  independent  variables  are  just  the  usual  kind  of  variables,  which  can  be  read  and 
/.agnod  to,  The  dofamdo'it  v,.iriat)les  are  virtual  variables  whose  values  are  defined  in  terms  of  the 
iiidepende''.t  van.Tjles.  rfio.se  values  can  be  read,  but  not  modified.  We  can  think  of  the  reading  of  a 
•:l';(i':ndont  va'iablo  as  shorthand  for  a  read  of  several  independent  variables,  together  with  a 
o.ilculation  of  the  lunction  giving  the  dependency. 

IndoponcJnnt  Vori.abics 

I'Evvl  ,'i  h  i  o,  f:;r  til  ’  Cut  o'  roguosts  that  originated  at  v, 

AC  t  'Vi'  f  ;,’  th,'  c‘'!  of  ’  .,;uu.  i'l  that  originated  at  v,  which  are  still  unsatisfied, 

;  n[  i;  ‘or  ..f ,  /  r  oor,  nc  'll  RESOURCES  that  have  not  yet  been  "captured"  by  requests, 

f  i'  oAf  i  r  ,r  1 .  f  u'.iiirod  ropource  on  its  way  back  to  a  request, 

APtUV.'C  ,v  I  Cf  If ,!  '■]<  OS.  for  the  number  of  "ARFtlVAL"  messages  received  from  each  of 

V  O  (  hi  hlr,.n  od  :  ot  to  v'c  'u  lit  rocpoc  llVu'ly, 
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s  €  bins,  then 

if  bin  s  is  nonempty, 

then  subtract  1  from  the  number  of  resources  in  s 
else  if  some  bin  is  nonempty  then 

[SELECT  the  first  unused  element  of  c  describing  a  nonempty  bin,  t; 
subtract  1  from  the  number  of  resources  in  t] 

X,  then 

if  some  bin  is  nonempty  then 

[SELECT  the  first  unused  element  of  c  describing  a  nonempty  bin,  t: 
subtract  1  from  the  number  of  resources  in  t] 

Define  SELECT(S,c.p,i)  to  be  the  number  of  times  bin  i  is  SELECTed  during  the  course  of 
processing  S  on  c  and  p.  (Note  that  a  bin  is  only  said  to  be  SELECTed  when  the  choice  sequence  is 
used  to  choose  it,  and  not  when  it  is  explicitly  chosen  by  the  script.  Define  choice(S,c,p,j)  to  be  equal 
to  k  provided  that  when  S(j)  is  processed  on  c  and  p,  the  kth  element  of  c  is  used  to  select  a  bin.  (If  no 
element  of  c  is  used,  then  choice(S,c,p.j)  is  undefined.)  It  follows  that  SELECT(S,c,p,i)  =  |{j: 
c(choice(S,c,p,j))  =  i}|, 

For  any  script  S,  let  binsequence(S)  denote  the  subsequence  of  S  consisting  of  bin  numbers. 
Script  S  is  said  to  dominate  script  S'  provided  that:  (a)  T  =  T’,  where  T  =  binsequence(S)  and  T'  = 
binsequence(S').  (b)  the  total  number  of  X's  in  S  is  at  least  as  great  as  the  total  number  o'  X's  in  S', 
and  (c)  for  each  i,  the  number  of  X's  in  S  preceding  T(i)  is  at  least  as  great  as  the  number  of  X's  in  S' 
preceding  T'(i).  The  main  result  of  this  section  is  that,  if  S  dominates  S',  then  SELECT(S,c,p,i)  > 
SELECT(S',c,p.i)  for  all  c,  p  and  i. 

We  say  that  an  interchange  of  two  consecutive  elements  of  a  script  S  is  legal  provided  that  the 
first  element  of  the  pair  is  an  X.  We  say  that  a  script  S’  is  reachable  from  a  script  S  if  S  can  be 
transformed  to  S'  by  a  series  of  legal  interchanges.  Note  that  S  dominates  all  scripts  S'  reachable 
from  Si  moreover,  if  S  dominates  S',  then  S'  can  be  augmented  with  some  suffix  of  X's,  to  a  script 
which  IS  reachable  from  S. 

Lemma  6:  For  any  scripts  S  and  S'  such  that  S'  is  reachable  from  S,  and  for  any 
clioice  seriuence  c,  [jiacement  function  p  and  bin  i, 

S[;L[:CT(S,c.p,i)  >  SELECT(S',c,p,i). 

Proof;  We  prove  this  lemma  by  showing  that  if  S'  is  reachable  from  S  by  a  single  legal 
interchange,  then  the  inequality  holds.  The  lemma  follows  by  induction  on  the  numlier  of 
leijal  interchanges. 

Fix  S,  c.  and  p.  Assume  that  S'  is  obtained  from  S  by  interchanging  S(j)  =  X  with 
:'.(i  I-  I)  If  S(|  t  I)  -  X.  then  S  S',  .so  the  result  is  obvious.  So  assume  S(i  +  1)  =  s  € 
bms  There  aor  thrp':'  (,.ises. 


C.K'O  1 :  Rin  s  is  omirty  after  processing  S(1)..  o(j  1)  on  c  and  p. 
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Lemma  4:  bnum-  =  vbnum-^. 

T,P  T,p 

Proof;  We  sketch  the  argument  for  fixed  r  and  C.  For  a  particular  edge  e.  let  a^  denote 

the  number  of  request  arrivals  below  e,  b^  denote  the  number  of  "BUYER"  messages  sent 

downward  along  e,  and  denote  the  number  of  resources  placed  below  e,  in  the 

execution  for  r  and  C.  Since  all  resources  get  matched  to  requests,  we  must  have  a^  +  b^ 

>  pg,  so  that  the  number  of  virtual  buyers  sent  over  edge  e  is  exactly  maxfa^  +  b^  -  p^,  0) 

=  a„  +  b„  •  p„. 
e  e 

Now  consider  all  the  edges  at  any  particular  height  h  in  'he  tree.  Since  all  resources 
and  requests  are  at  the  leaves,  and  the  branches  are  all  of  equal  length,  it  is  clear  that 

^heigh:.^(e)  =  h  ^  total(p)  =  ^height^{e)  =  h  ^e' 

Therefore, 

^heiyiit^(e)  -  h  “height.j.(e)  =  h  ^ 

Thct  s,  the  numbers  of  buyers  and  virtual  buyers  sent  over  edges  of  height  h  are  equal. 
Since  this  IS  true  for  all  h,  the  result  follows. I 
Theorem  5;  COStj.^  <  4(bnum.|.p). 


Proof:  Iminodiato  from  Lemmas  3  and  Restriction  3.1 


rtius  iM  Ollier  to  obtain  an  upper  bound  on  costj.  ^ff),  it  suffices  to  prove  a  bound  for  bnum.|.  ^(1). 

4  4  A  Combinatorial  Result 

I  ii'S  CLib'Vctinn  contains  a  key  combinatorial  result  which  will  be  used  in  the  subsequent  analysis, 
'ni  f;..’!  :  navior  of  the  algorithm  at  a  single  node  V.  The  children  of  v  are  modelled  as  a  set  of 

;  ’ .  for  resoLirces.  (Here,  we  rio  not  concern  ourselves  about  the  tree  structure  beyond  the  children.) 
i  •  t  (,  be  a  ch,  ice  ;:,oquence  for  v  Erich  bin  s  is  initialized  to  contain  a  number  p(s)  of  resources. 

I  he  arrival  of  messages  at  v  is  described  by  a  script,  S.  A  script  is  a  finite  sequence  of  symbols, 
■ir.fi  of  v.hicfi  '$  either  a  bin  number  s  or  an  "X".  A  bin  number  represents  the  arrival  of  an 
'AHb'lVAL  '  ir’'r,r,age  from  the  specified  child.  The  symbol  X  represents  the  arrival  of  a  "BUYER" 
nu.  s  ',.'ge  from  v's  parent. 

hli.,'  prooe  .Miig  of  script  S  on  c  and  p,  is  as  follows.  The  elements  of  S  are  processed 

:;i)'j*;ntially. 

if  ;■> (  1  )  is: 


Lemma  3:  (a'  searchcost^  <  bnum^  +  vbnum,.  . 

1  .p  ““  I  ,p  T.p 

(b)  returncost,.  ^  =  gnum-j^  p. 

(c)  gnum^  p  <  bnum^  p  +  vbnum^  p. 

(d) cost^p  <  2[bnum^p  +  vbnum^p]. 

Proof:  (a)  This  inequality  is  true  because  buyers  continue  to  make  progress  up  and 
down  the  edges  of  the  tree;  all  time  used  by  the  algorithm  is  occupied  by  the  transmission 
of  appropriate  buyer  and  virtual  buyer  messages.  The  reason  that  we  have  an  inequality 
rather  than  an  equation  here  is  that  buyers  are  permitted  to  travel  "discontinuously",  as 
described  in  Section  3. 

(b)  This  equation  is  true  because  captured  resources  travel  continuously  via  captured 
resource  messages. 

(c)  We  must  show  that  each  captured  resource  message  always  moves  in  such  a  way 
as  to  "negate"  a  buyer  or  virtual  buyer  message.  This  is  a  bit  tricky  to  argue,  because  of 
the  discrepancies  between  estimates  at  opposite  ends  of  an  edge. 

A  captured  resource  only  moves  over  an  edge  if  the  net  flow  of  buyers  into  the  node  on 
that  edge,  as  estimated  at  the  near  endpoint,  is  positive.  By  moving  over  that  edge,  the 
captured  resource  negates  an  incoming  buyer  or  virtual  buyer  along  that  edge,  as 
estimated  at  the  near  endpoint.  Because  of  the  assumption  that  all  messages  fake  exactly 
time  1,  by  the  time  the  captured  resource  reaches  the  far  endpoint,  the  negated  buyer  or 
virtual  buyer  is  also  counted  in  the  estimate  of  outgoing  buyers  and  virtual  buyers  at  the 
opposite  endpoint.  The  arrival  of  the  captured  resource  at  the  far  endpoint  can  thus  be 
regarded  as  negating  an  outgoing  buyer  or  virtual  buyer  at  the  far  endpoint  as  well. 

(d)  Straightforward  by  Lemma  2  and  (a)  (c).l 

Now,  we  introduce  an  additional  restriction,  to  remain  in  force  for  the  remainder  of  Section  4. 

Restriction  3: 

T  has  all  leaves  at  the  same  distance  from  the  roof,  and  r  and  p  are  nonzero  only  at  leaves. 

As  usual,  the  following  lemma  is  intended  to  hold  for  all  valid  domains  of  definition. 
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Let  f  denote  an  arbitrary  probability  density  function  whose  domain  consists  of  positive  reals.  Extend 
the  domain  of  the  function  cost.^  p  to  the  set  of  such  functions  by  defining  costj^lf)  to  be  the  expected 
value  of  cost.|.p(r),  where  r  is  of  length  total(p),  with  its  successive  locations  chosen  independently 
using  the  distribution  <pj,  and  its  successive  interarrival  times  chosen  using  f.  That  is,  at  the  time  the 
algorithm  begins,  and  at  the  time  of  each  request,  the  probability  that  the  next  arrival  occurs  exactly  t 
units  later  is  f(t).  We  will  be  primarily  interested  in  this  cost,  cost.j.  p(f). 

We  define  searchcostj  ^(r),  searchcostj  ^(1) ,  etc.,  analogously  to  our  earlier  definitions. 

The  following  claim  is  true  for  all  domains  for  which  the  definitions  are  valid. 

Lemma  2;  cost,  „  =  searchcost,  „  +  returncost, 

T.P  T.p  T,p 

Proof:  Straightforward.! 

Next,  we  will  relate  the  given  cost  measures  to  the  total  numbers  of  various  kinds  of  messages  sent 
during  the  execution  of  the  algorithm.  Note  that  during  the  execution  of  an  algorithm,  the  estimates 
of  "BUYER"  and  virtual  buyer  messages  sent  along  an  edge  can  be  different  at  the  two  ends  of  the 
edge.  However,  after  the  entire  execution  of  the  algorithm  is  completed,  the  discrepancy  disappears, 
so  that  the  following  definitions  are  unambiguous.  Let  bnumj^(r,C)  denote  the  total  number  of 
"BUYER"  messages  sent  on  all  edges  during  the  execution  of  the  algorithm  on  r  using  C.  Let 
vbnurrij  Jr,C)  denote  the  total  number  of  virtual  buyer  messages  sent  on  all  edges  during  the 
execution  of  the  algorithm  on  r  using  C.  Let  gnunij  J(r,Z)  denote  the  total  number  of  captured 
resource  messages  sent  on  all  edges  during  the  execution  of  the  algorithm  on  r  using  C.  As  before, 
define  bnunij. ^(r),  bnum^^(l) ,  etc. 

Because  of  the  fact  that  message  delivery  time  is  assumed  to  be  exactly  1,  there  are  some 
relationships  between  the  measures  describing  time  costs  and  the  measures  describing  numbers  of 
messages.  The  following  lemma  describes  a  set  of  relationships  among  the  various  measures.  Note 
that  all  the  statements  are  true  over  all  possible  domains  of  definition. 
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Restriction  2  has  the  effect  of  restricting  the  executions  under  consideration;  for  eornple,  all 
messages  between  any  two  nodes  are  pipelined  -  they  arrive  in  the  order  in  which  they  are  sent. 
While  we  would  like  to  understand  the  behavior  of  the  algorithm  in  the  presence  of  variable  message 
delivery  times,  such  analysis  appears  to  be  more  difficult. 

4.3.  Cost  Measures  and  Preliminary  Results 

A  request  pattern,  r,  is  a  finite  sequence  of  elements  of  vertices.|.  x  R  whose  second  components 
are  monotone  nondecreasing.  A  request  pattern  represents  the  sequence  of  requests  that  occur, 
their  locations  and  times.  If  r  is  a  request  sequence,  then  length(r)  denotes  its  length. 

A  choice  sequence,  c,  for  v  €  internal.j.  is  an  infinite  sequence  of  elements  of  childron.j.(v),  with 
infinitely  many  occurrences  of  each  child.  If  C  =  {c^}  is  a  collection  of  choice  sequences,  one  for 
each  V  €  internal.^,  then  C  can  be  used  in  place  of  probabilistic  choices  in  an  execution  of  the 
algorithm,  as  follows.  Each  internal  node,  v,  makes  choices  among  its  children  by  choosing  the  first 
unused  element  of  c^  satisfying  the  inequality  PLACE(desc^(s))  >  ARRIVALS(s)  +  BUYERS(s).  That 
is,  V  chooses  a  child,  s,  for  which  v  thinks  there  are  still  remaining  resources  in  s’s  subtree. 

Let  p  be  a  placement  for  T.  Let  r  be  a  request  pattern,  and  C  =  {c^}  a  collection  of  choice 
sequences,  one  for  each  v  €  internal.|..  Then  costjjir.t)  is  defined  to  be  the  total  time  from  requests 
to  corresponding  grants,  if  requests  arrive  according  to  r  and  C  is  used  in  place  of  probabilistic 
choices.  (With  suitable  conventions  for  handling  events  which  happen  at  the  same  time,  the 
execution,  and  hence  the  cost,  is  uniquely  defined  for  fixed  r  and  C.) 

The  cost  measure  defined  above  can  be  broken  up  into  two  pieces,  as  follows.  Let 
scarchcost^  ^(r,0  be  the  total  of  the  times  from  requests  to  corresponding  captures  of  resources,  if  r 
and  C  are  used  as  above.  Let  rcturncostj ^(r,C)  be  the  total  of  the  times  from  captures  to 
corresponding  grants  of  resources. 

Now  we  incorporate  a  probabilistic  construction  of  C  into  the  cost  measure.  If  r  is  a  request 
pattern,  lot  cost ^  Jr)  denote  the  expected  value  of  cost.|.  ^^(r,C),  whore  C  is  constructed  using  (pj.. 
(That  is,  for  each  v  C  intornal|,  the  sequence  c^  is  constructed  by  successive  choices  from  among 
children^(v),  choosing  s  v;ilh  probability  (desc., (s))/(f)^(desCj^(S)),  where  S  =  childrcn^(v).  Among 
the  sequences  thereby  generated  are  some  for  which  it  is  not  the  case  that  each  child  occurs 
infinitely  often.  However,  these  sequences  form  a  set  of  measure  0,  so  that  we  can  ignore  them  in 
calculating  the  expected  cost  measure.)  We  claim  that  cost.|.p(r)  is  exactly  the  expected  total  time 
fiom  reqiinsts  to  (jr.ints,  providi.’d  thr;  algorithm  is  run  in  the  normal  way,  using  prnb.abili'.tic  choices. 
That  is,  the  two  strategies  of  constructing  choice  sequences  independently  of  the  algorithm  and 
carrying  out  the  probabilistic  choices  on  line  give  identical  results. 
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If  V  is  a  vertex  of  T,  let  height  jiv)  denote  the  maximum  distance  from  v  to  a  leaf  in  its  subtree.  If  e  is  an 
edge  in  T,  then  define  height j(e)  to  be  the  same  as  height.|.(v),  where  v  is  e’s  upper  endpoint.  Let 
height j  denote  height.j.(root.j.). 

A  placement  for  T  is  a  function  p:  vertices^  — »  N,  representing  the  number  of  resources  at  each 
vertex.  We  write  total(p)  tor  p(vertices.j.),  the  total  number  of  resources  in  the  entire  tree.  We  say  that 
p  is  nonnull  provided  total(p)  >  0. 

A  weighted  tree,  T,  is  an  undirected,  rooted  tree  with  an  associated  probability  density  function, 

on  the  leaves  of  T,  such  that  <p^(v)  >  0  for  all  leaves  v.  (This  assumption  is  made  for  technical 
reasons,  so  that  we  can  normalize  probability  functions  without  danger  of  dividing  by  0.)  if  T  is  a 
weighted  tree,  v  €  internal.^,  and  S  is  a  nonempty  subset  of  children.|.(v),  then  let  random  j  ^  denote  the 
probability  function  which  returns  s  €  S  with  probability  <p.|.{desc.|.(s))/<p.|.(desc.^(S)).  Thus,  random.j.g 
returns  s  with  probability  proportional  to  the  sum  of  the  probability  function  values  for  the 
descendants  of  s. 

4.2.  Initial  Restrictions 

For  the  remainder  of  Section  4,  we  assume  that  the  following  two  restrictions  hold. 

Restriction  1: 

T  is  a  weighted  tree,  and  the  nondeterministic  choice  step  in  Part  (2)  of  the  algorithm  uses  a  call  to 
random.|.g, 

Restriction  2: 

Delivery  time  for  messages  is  always  exactly  1. 

Restriction  1  describes  a  particular  method  of  choosing  among  alternative  search  directions.  This 
rnettiod  does  not  use  all  the  information  available  during  execution,  but  only  the  "static"  probability 
distribution  information  available  at  the  beginning  of  execution.  One  might  expect  a  more  adaptive 
choice  method  to  work  better;  however,  we  do  not  know  how  to  analyze  such  strategies. 
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be  the  number  of  captured  resource  messages  which  have  been  sent  from  w  to  v  but  have 
not  yet  arrived.  (Of  course,  neither  v  nor  w  actually  "knows"  the  value  of  this  vari.ible.) 
For  any  time,  t,  after  the  NETBUYERS(w)  and  REQUEST  values  have  stabilized,  any  node 
V,  and  any  w  €  NEIGHBORS(v),  let  A(v,w,t)  be  the  value  of  NETBUYERS  w)  - 
NETGRANTS(w)  +  MESSAGES(w)  at  v  at  time  t.  Note  that  A(v,w,t)  =  -Afw.v.t)  in  all  cases. 
Let  SUM(t)  denote  ^|A(v,w,t)|.  We  claim  that  any  event  which  involves  the  receipt  of  a 
captured  resource  message  does  not  change  SUM(t),  while  any  event  which  involves  the 
sending  of  a  captured  resource  message  decreases  SUM(t).  Therefore,  captured  resource 
messages  will  not  be  sent  forever:  they  will  eventually  subside,  at  which  time  they  must 
have  found  a  matching  request. 

First,  consider  an  event  involving  the  receipt  of  a  captured  resource,  by  v,  from  w.  The 
only  term  in  the  sum  which  is  affected  is  A(v,w,t).  The  receipt  of  the  messages  causes  v’s 
values  of  MESSAGES(w)  and  NETGRANTS(w)  both  to  decrease  by  1,  so  that  A(v,v/,t)  is 
unchanged.  Therefore,  SUM(t)  is  unchanged.  Second,  consider  an  event  involving  the 
sending  of  a  captured  resource,  by  v,  to  w.  The  only  terms  in  the  sum  which  are  affected 
are  A(v,v/.t)  and  A(w,v,t).  At  time  t  just  prior  to  the  sending  event,  it  must  be  that  v's  value 
of  NETBUYERS(w)  -  NETGRANTS(w)  >  0,  which  implies  that  A(v,w,t)  >  0.  The  result  of 
sending  the  message  is  to  increase  NETGRANTS(w),  which  means  that  A(v,w,t)  gets 
decreased  by  1.  Therefore,  |A(v,w,t)l  gets  decreased  by  1.  Thus,  also,  |A(w,v,t)|  gets 
decreased  by  1 ,  so  that  SUM(t)  gets  decreased  by  2.1 


4.  Monotonicity  Analysis 

The  rest  of  the  paper  is  devoted  to  an  analysis  of  the  time  requirements  of  the  algorithm. 
Specifically,  we  measure  the  sum  of  the  times  between  requests  and  their  corresponding  g.’-ants.  For 
the  purpose  of  carrying  out  the  analysis,  certain  restrictions  will  be  made.  These  restrictions  will  be 
introduced  as  needed. 

We  begin  with  some  basic  definitions.  Next,  we  introduce  two  restrictions  which  are  needed 
throughout  the  analysis.  Then  we  define  and  categorize  the  complexity  measures  of  interest.  We 
then  prove  a  basic  combinatorial  result,  and  use  it  to  prove  the  monotonicity  of  the  number  of 
"DUYEIT'  messages  as  a  function  of  interarrival  time.  Finally,  we  show  that  the  expected  running 
time  of  the  .algorithm  is  bounded  by  the  expected  time  for  the  sequential  case  of  the  algorithm. 


4. 1 .  Definitions 

Lot  N  denote  the  sot  of  natural  numbers,  including  0.  Let  R"^  denote  the  set  of  nonnegtitive  reals. 
If  f  is  a  numerical  function  with  domain  V,  then  extend  f  to  subsets  of  V  by  f(W)  = 

Let  T  be  a  rooted  tree.  We  write  vertices j.,  intcrnnij,  and  IsnveSj  to  denote  the  indicated  sets  of 
vertices  of  T.  Lot  rootj  denote  the  root.  If  v  G  verticeSp  we  write  desCj(v)  for  the  set  of  vertices  of  T 
which  are  descendants  of  v  (including  v  itself),  parent for  v's  parent  in  T,  children j(v)  for  v’s 
children,  and  neicjhhoir, for  childron^(v)  U  (parent |.(v)}. 


This  expression  is,  in  turn,  equal  to 


PLACE(DESCENDANTS)  ■  5:^£j,^,^qp^^PLACE(DESCENDANTS(s)), 

=  PLACE(v)  =  |CAPTURED|. 

Thus,  NETFLOW  <  |CAPTURED|,  a  contradiction. 

Thus,  we  have  shown  that  it  is  always  possible  to  service  an  excess  request. 

Next,  we  must  show  that  NETFLOW  =  |CAPTURED)  between  Parts  2  and  3  of  the  code. 
This  means  that  after  servicing  any  excess  request,  there  is  no  remaining  request  to  be 
serviced.  Previous  to  Part  2,  |CAPTURED|  <  NETFLOW  <  |CAPTUREDl  +1.  If 
NETFLOW  was  equal  to  |CAPTUREDl  +  1,  then  the  body  of  the  conditional  was  executed. 
If  the  first  case  of  the  conditional  held  (i.e.  the  case  for  FREE  ^  0),  then  |CAPTURED|  is 
increased  by  1,  so  the  invariant  is  restored.  Otherwise,  a  "BUYER"  message  was  sent  to  a 
child,  s,  for  which  PLACE(DESCENDANTS(s))  -  ARRIVALS(s)  >  BUYERS(s).  This  caused 
NETBUYERS(s)  to  increase  by  t ,  thereby  increasing  the  value  of  NETFLOW  and  restoring 
the  invariant. 

The  third  portion  of  the  algorithm  manages  the  travel  of  captured  resources  back  to 
requests  First,  note  that  there  can  be  only  one  captured  resource  assigned  to  GRANT  at 
any  node  in  a  single  step,  since  the  two  assignments  to  GRANT  cannot  both  be  executed 
during  a  single  step.  If  the  message  is  a  captured  resource,  then  no  progress  is  done  until 
tfie  clause  contains  the  second  grant.  Otherwise,  this  clause  is  skipped.  We  must  argue 
that  such  a  neighbor  exists  in  this  case. 

Assume  not.  Then  NETBUYERS(s)  <  NETGRANTS(s)  for  all  s  €  NEIGHBORS.  Now, 
NETFLOW  =  ICAPTUREDI,  so  that  ICAPTURED)  =  |REQUESTS|  +  NETBUYERS, 

<  |REQUESTS|  +  NETGRANTS, 

=  IREQUESTSl  +  ICAPTURED]  -  jSATlSFlEDj  •  1 

=  lACTIVE]  +  ICAPTURED]  -  1.  Therefore,!  <  jACTIVEl,  a  contradiction. 

Thus,  we  have  checked  that  the  key  assertions  hold  and  the  code  can  be  executed  at 
ail  points.  We  have  claimed  (and  tried  to  argue)  that  the  algorithm  follows  the  strategy  of 
the  preceding  section,  in  setting  up  a  flow  of  buyers  from  requests  to  resources. 
Eventually,  the  values  of  all  the  NETBUYERS(w)  variables  will  stabilize,  and  the  values 
taken  on  by  corresponding  NETBUYERS(w)  variables  at  either  end  of  a  single  edge  will  be 
negations  of  each  other.  (We  use  the  fact  that  there  are  only  finitely  many  requests  here. 
Eventually,  no  further  requests  will  arrive,  so  no  additional  "ARRIVAL"  messages  will  be 
sent.  There  is  a  bound  on  how  many  "BUYER"  messages  will  be  sent  dov/nward  along 
any  edge.  Therefore,  there  are  only  finitely  many  total  "ARRIVAI."  and  "BUYER" 
messages  which  get  sent,  so  that  eventually,  they  will  all  be  delivered.)  Similarly,  all  the 
REQUESTS  variables  will  eventually  stabilize. 

Finally,  we  must  consider  the  travel  of  captured  resources  to  request  origins.  Define  a 
new  variable,  MESSAGESfw),  at  node  v,  where  w  G  NEIGHBORS(v).  Its  value  is  defined  to 
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By  the  ARRIVAL  invariant,  this  is  equal  to 

PLACE(DESCENDANTS(s))  -  I.^children,  t^s 


However,  the  original  invariant  says  that  NETFLOW  = 
is  never  permitted  to  be  greater  than  PLACE(v),  a 

We  have  thus  shown  that  |CAPTURED|  <  NETFLOW  <  |CAPTURED|  +  1  at  the  point 
where  that  claim  is  made.  Thus,  there  is  at  most  one  excess  request  that  requires 
disposition.  In  the  case  where  there  is  an  excess  request,  node  v  must  service  that  request 
in  its  subtree.  There  are  two  possibilities:  either  v  can  service  the  request  locally,  or  it 
cannot.  If  FREE  0,  then  a  free  resource  is  captured  to  service  the  excess  request.  If 
not,  then  a  "BUYER"  message  must  be  sent  down  into  some  subtree.  We  must  show  that, 
in  the  event  FREE  =  0,  it  is  possible  to  send  such  a  "BUYER"  message.  That  is,  we  must 
check  that  S  0  at  the  place  where  that  claim  is  made. 

Assume  not.  We  will  make  some  deductions  about  the  values  of  the  variables  at  the 
point  where  that  claim  is  made.  At  that  point,  we  know  that  FREE  =  0,  so  that  PLACE(v) 
=  |CAPTURED|.  We  also  know  that 

NETBUYERS(PARENT)  <  PLACE(DESCENDANTS)  •  ARRIVALS(PARENT),  by 
definition  of  NETBUYERS. 

Then 

NETFLOW  =  |REQUESTS|  +  NETBUYERS(PARENT)  +  5;^^^^,^Pf,g^NETBUYERS(s), 

<  IREQUESTSI  +  PLACE(DESCENDANTS)  ■  ARRIVALS(PARENT)  + 

Because  S  =  0,  it  follows  that 

PLACE(DESCENDANTS(s))  -  ARRIVALS(s)  <  BUYERS(s)  for  each  s  E  CHILDREN. 

Therefore, 

NETBUYERS(s)  =  (PLACEfDESCENDANTSis))  ■  ARRIVALS(s)). 

Thus,  the  right-hand  side  of  the  next-to-last  inequality  is  equal  to 

IREOUESTSI  +  PLACE(DESCENDANTS)  ARRIVALS(PARENT) 

^'sEchildren(P'-^^^(DESCENDANTS(s))  -  ARniVALS(s)).  - 


PLACE(DESCENDANTS)  • 
PLACE(DESCENDANTS(t)), 

=  PLACE(v). 

Thus.  NETFLOW  >  PLACE(v). 
ICAPTUREDI,  and  |CAPTURED| 
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The  second  portion  of  the  algorithm  manages  the  disposition  of  any  excess  flow  of 
requests  into  the  node.  We  must  first  check  that  the  number  of  excess  requests  after  the 
initial  processing  of  a  single  message  can  only  be  0  or  1.  That  is,  we  must  verify  that 
|CAPTURED|  <  NETFLOW  <  |CAPTURED|  +  1  between  Parts  1  and  2  of  the  code. 

A  quick  check  of  the  cases  shows  that  the  only  way  this  could  fail  to  be  true  is  if  M  = 
"ARRIVAL"  from  a  child  s,  and  the  result  of  processing  M  causes  NETBUYERS(s)  to 
remain  unchanged,  while  NETBUYERS(PARENT)  decreases.  In  this  case,  we  can  deduce 
some  relationships  among  the  values  of  v’s  local  variables  at  the  beginning  of  the  node 
step. 

For  every  t  €  CHILDREN,  it  must  be  the  case  just  before  execution  of  Part  1  that 
-NETBUYERSft)  =  min(PLACE(DESCENDANTS(t))  -  ARRIVALS(t),  BUYERS(t)), 
so  that 

-NETBUYERS(t)  <  PLACE(DESCENDANTS(t))  ■  ARRIVALS(t). 

That  is, 

NETBUYERS(t)  >  ARRIVALS(t)  ■  PLACE(DESCENDANTS(t)). 

Since  NETBLIYERS(s)  remains  unchanged,  then  it  must  be  the  case  that 

-NETBUYERS(s)  PLACE(DESCENDANTS(s))  •  ARRIVALS(s). 

(If  they  were  equal,  then  PLACE(DESCENDANTS)(s))  •  ARRIVALS(s)  * 
min(PLACE(DESCENDANTS(s))  -  ARRIVALS(s),  BUYERS(s)),  and  an  increase  to 
AnRIVALS(s)  would  cause  a  change  to  the  minimum,  thereby  changing  NETBUYERS(s).) 

Therefore,  NETBUYERS(s)  ARRIVALS{s)  •  PLACE(DESCENDANTS(s)),  and  so 

NETBUYERS(s)  >  ARRIVALS(s)  ■  PLACE(DESCENDANTS(s)). 

Since  NETBUYERS{PARENT)  decreases,  it  means  that  NETBUYERS(PARENT)  = 
PLACE(DESCENDANTS)  -  ARRIVALS(PARENT). 

Now  consider  NETFLOW  =  |REQUESTS|  +  NETBUYERS.  The  right  side  is  equal  to 

IREQUESTSI  +  NETBUYERS(PARENT)  +  NETBUYERS(s)  + 
,^gNETBUYERS(w). 

By  previous  results,  this  is,  in  turn,  strictly  greater  than 
IREQUESTSI  +  PLACE(DESCENDANTS)  ■  ARRIVALSfPARENT) 

+  ARRIVALS(s)  -  PLACE(DESCENDANTS(s)) 

^  ^•teCMiLDHi  M  tv^sARRIVALStt)  ■  PLACE(DESCENDANTS(t)). 


(Part  3) 

/*  Process  M  if  M  is  a  captured  resource  message.  •/ 

If  M  is  a  resource  from  w  then 

[NETGRANTS{w)  :=  NETGRANTS(w)  -  1 
GRANT  :=  M] 

/•  Send  a  captured  resource,  if  you  have  one,  toward  a  request  origin.  •/ 

If  GRANT  ^  0  then 

/*  NETGRANTS  =  |CAPTURED|  -  ISATISFIEOj  -  1.  •/ 

[if  ACTIVE  ^  0  then 
then 

[choose  r  €  ACTIVE 
ACTIVE  :=  ACTIVE  -  (r) 
output  (r, GRANT)] 
else 

[choose  w  €  NEIGHBORS  with  NETBUYERS(w)  -  NETGRANTS{w)  >  0 
send  GRANT  to  w 

NETGRANTS(w)  ;=  NETGRANTS(w)  +  1] 

GRANT  :=  0] 

3.  Correctness  of  the  Algorithm 

Theorem  1 :  The  given  algorithm  solves  the  resource  allocation  problem. 

Proof:  Vi/e  claim  that  the  node  program  given  above  implements  the  strategy 
described  informally  in  the  previous  section.  We  do  not  give  a  proof  of  this 
correspondence  here.  Rather,  we  argue  correctness  of  the  key  assertions  of  the  program 
and  give  informal  arguments  for  the  rest  of  the  proof  of  correctness  of  the  algorithm. 

The  first  portion  of  the  algorithm,  the  initial  processing  of  the  first  three  kinds  of 
messages,  simply  sends  the  appropriate  "ARRIVAL"  messages  and  records  the  proper 
changes  to  the  various  sets  and  counters. 

For  any  of  the  three  kinds  of  messages,  node  v  is  finding  out  about  a  new  request  that 
needs  to  be  processed.  In  some  cases,  v  will  need  to  do  more  to  help  process  the  request. 
If  the  message  is  an  "ARRIVAL",  and  node  v  thinks  that  the  corresponding  request  can  be 
serviced  in  the  sender's  subtree,  then  v  has  no  further  work  to  do.  If  the  message  is  a 
request  or  an  "ARRIVAL",  and  if  node  v  thinks  that  it  is  impossible  to  service  that  request 
in  v's  subtree,  then  the  "ARRIVAL"  message  sent  upward  by  v  will  be  counted  by  v’s 
parent  in  its  estimate  of  virtual  buyers  emanating  from  v’s  subtree.  Thus,  after  sending  this 
"ARRIVAL"  message  upward,  v  will  have  no  further  work  to  do.  Also,  if  the  message  is  a 
"BUYER"  and  v  thinks  that  it  is  impossible  to  service  the  request  in  v’s  subtree,  then  v 
need  not  do  anything  more.  However,  if  v  thinks  that  the  new  request  can  be  serviced  in  its 
subtree,  then  it  has  some  further  work  to  do,  in  the  second  portion  of  the  algorithm. 
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/•  The  following  invariants  hold  at  the  beginning  of  any  node  step. 

NETFLOW  =  ICAPTUREDl . 

ARRIVALS(PARENT)  =  IREQUESTSl  5;^  £  (.„jldr£„ARRIVALS( w) ] . 

GRANT  =  0. 

NETGRANTS  =  ICAPTUREOj  -  jSATISFIEOl.  •/ 

{Part  1) 

If  M  is  a  request  then 

[REQUESTS  :=  REQUESTS  U  {M} 

ACTIVE  :=  ACTIVE  U  {M} 
send  "ARRIVAL"  to  PARENT 
ARRIVALS(PARENT)  :=  ARRIVALS( PARENT)  +  1] 

If  M  =  "ARRIVAL"  from  w  then 
[ARRIVALS(w)  ;=  ARRIVALS(w)  +  1 
send  "ARRIVAL"  to  PARENT 
ARRIVALS(PARENT)  :=  ARRIVALS(PARENT)  +  1] 

If  M  =  "BUYER"  then  BUYERS( PARENT)  :=  BUYERS( PARENT)  +  1 

/*  Now  slightly  revised  invariants  hold: 

ICAPTUREDl  <  NETFLOW  <  |CAPTURED|  +  1. 

ARRIVALS(PARENT)  =  IREQUESTSl  +  \  £  j.„jldrjhARRIVALS(w) ] . 

GRANT  =  0. 

NETGRANTS  =  |CAPTUREO|  -  |SATISFIE01.  */ 

(Part  2) 

/*  Next,  if  there  is  an  excess  request,  service  it.  */ 

If  NETFLOW  =  ICAPTUREDl  1  then 

if  FREE  ^  0 
then 

[choose  s  €  FREE 
FREE  ;=  FREE  -  {s} 

GRANT  :=  Ss] 
el  se 

[/*  Send  "BUYER"  down  into  a  subtree.  ♦/ 

S  :=  (s  e  CHILDREN:  PLACE(DESCENDANTS(s) )  >  ARRIVALS(s)  +  BUYERS(s)} 

/*  S  ^  0.  •/ 

choose  NEXT  €  S 
send  "BUYER"  to  NEXT 
BUYERS(NEXT)  ; BUYI:RS( NEXT )  +  1] 


/*  NETFLOW  =  |CAl>rURCD|.  */ 
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Another  way  to  understand  the  equations  is  as  follows.  Again,  consider  the  first  equation,  for  w  = 
PARENT.  Then  NETBUYERS(w)  =  BUYERS(w)  ■  VIRTBUYERS(w),  where  the  latter  quantity  is  the 
number  of  virtual  buyers  which  v  estimates  it  has  sent  to  its  parent.  Using  the  expression  which  was 
derived  in  the  preceding  subsection  for  the  number  of  virtual  buyers,  we  see  that  NETBUYERS(w)  = 
BUYERS(w)  ■  max(ARRIVALS(w)  +  BUYERS(w)  -  PLACE(DESCENDANTS),0).  This  is  equal  to 
min(PLACE(DESCENDANTS)  -  ARRIVALS(w),BUYERS(v.  „  as  needed.  Again,  the  other  calculation 
is  similar. 

The  remaining  dependent  variables  are: 

NETBUYERS,  for  the  total  of  all  the  NETBUYERS(w), 

Dependency:  NETBUYERS  =  51^  £  i^^iQ^gQpjgNETBUYERSfw). 

NETFLOW,  for  the  net  flow  of  buyers  into  v. 

Dependency:  NETFLOW  =  |REQUESTS|  +  NETBUYERS. 

NETGRANTS,  for  the  net  flow  of  grants  out  of  v, 

Dependency:  NETGRANTS  =  £  |,^giQ^gQP,gNETGRANTS(w). 

The  following  code  is  executed  in  response  to  the  receipt  of  any  message  or  request,  M.  The  first 
part  of  the  code  does  initial  processing  of  messages,  updating  estimates  and  sending  any  required 
"ARRIVAL"  messages.  After  the  first  part  of  the  code,  there  will  be  at  most  one  excess  request  left  at 
the  node,  and  if  there  is  such  a  request  left  at  the  node,  then  the  node  is  able  to  service  that  request, 
either  locally  or  by  sending  a  buyer  into  a  subtree.  In  the  second  part  of  the  code,  the  node  decides 
where  it  can  service  such  an  excess  request,  and  it  does  so.  (In  case  a  buyer  is  sent  down  into  a 
subtree,  the  subtree  is  chosen  nondeterministically.  Later,  we  will  refine  the  algorithm  to  use  a 
probabilistic  choice  at  this  point.)  Finally,  in  the  third  part  of  the  code,  the  node  processes  an  excess 
captured  resource,  if  it  happens  to  have  one.  (It  cannot  have  more  than  one.)  The  node  can  have  a 
captured  resource  either  because  M  was  a  captured  resource  message,  or  because  (in  the  second 
part  of  the  code)  the  node  itself  captured  a  local  resource.  The  resource  is  granted  to  a  local  request 
if  possible;  otherwise,  it  is  sent  in  such  a  way  as  to  negate  the  net  flow  of  buyers  into  the  node  along 
some  edge.  It  is  always  possible  to  process  such  a  captured  resource  in  one  of  these  ways. 
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BUYERS(w),  w  £  NEIGHBORS,  for  the  number  of  "BUYER"  messages  sent  to  v's  children  and 
received  from  v's  parent,  respectively, 

NETGRANTS(w),  w  £  NEIGHBORS,  for  the  net  flow  of  captured  resources  out  of  v  to  each  of  its 
neighbors, 

NEXT,  a  temporary  variable  which  can  hold  a  vertex. 

Initialization  of  Independent  Variables 

REQUESTS  =  ACTIVE  =  0,  FREE  =  RESOURCES,  and  all  other  variables  are  0. 

Dependent  Variables  and  their  Dependencies 

CAPTURED,  for  the  set  of  resources  in  RESOURCES  which  have  been  captured, 

Dependency;  CAPTURED  =  RESOURCES  -  FREE. 

SATISFIED,  for  the  set  of  requests  in  REQUESTS  which  have  been  satisfied, 

Dependency:  SATISFIED  =  REQUESTS  -  ACTIVE. 

NETBUYERS(w),  w  £  NEIGHBORS,  for  the  net  flow  of  buyers  and  virtual  buyers  into  v  from  each 
neighbor,  (Recall  the  definition  of  "virtual  buyers"  from  the  last  subsection.) 

Dependency:  If  w  =  PARENT,  then  NETBUYERS(w)  =  min(PLACE(DESCENDANTS)  ■ 
ARRIVALS(w),  BUYERS(vv)),  If  w  £  CHILDREN,  then  NETBUYERS(w) 
•min(PLACE(DESCENDANTS(w)  -  ARRIVALS(w).  BUYERS(w)). 

These  two  equations  can  be  understood  as  follows.  Consider,  for  example,  the  first  equation,  for 
vv  =  PARENT.  If  PLACE(DESCENDANTS)  -  ARRIVALS(w)  <  BUYERS(w),  it  means  that  the  placement 
originally  given  for  v's  subtree  is  not  adequate  for  handling  the  requests  (arrivals)  which  have 
originated  in  v's  subtree,  together  with  the  "BUYER"  messages  sent  down  from  w.  Therefore,  all  the 
resources  in  v's  subtree  are  allocated  to  requests,  either  within  or  outside  of  v's  subtree.  Whether  the 
net  flow  of  buyers  should  be  regarded  as  into  or  outward  from  v's  subtree  then  depends  solely  on  the 
.sign  of  PLACE(DESCENDANTS)  •  ARniVALS(w),  without  regard  to  the  number  of  "BUYER" 
messages  received  from  w.  That  is,  if  PLACE(DESCENDANTS)  <  ARRIVALS(w),  then  the  sign  is 
negative  and  the  net  flow  of  buyers  is  outward  from  v’s  subtree,  while  otherwise  it  is  inward;  in  either 
case,  its  magnitude  is  equal  to  |PLACE(DESCENDANTS)  •  ARRIVALS(w)|.  On  the  other  hand,  if 
PLACE(Dl;SCENDANTS)  ■  ARRIVALS(w)  >  BUYERS(w),  then  the  placement  originally  given  for  v's 
subtree  is  adequate  for  handling  both  the  requests  which  have  originated  in  v’s  subtree,  together  with 
the  "BUYER"  messages  sent  down  from  w.  Therefore,  the  net  flow  of  buyers  is  inward,  and  its 
amount  is  just  equal  to  BUYEnS(w),  v/ithout  regard  to  the  other  two  values.  The  second  equation  is 
siinil,ir,  with  a(ipropriate  changes  of  sign. 


Then  choices  using  c  are  made  for  both  S  and  S'  at  both  steps  j  and  j  +  1.  Thus, 
choice(S,c,p,j)  =  choice(S',c,p,j)  and  choice(S,c,p,j  + 1)  =  choice(S’,c,p,j  + 1).  The 
number  of  resources  remaining  in  each  bin  after  step  j  +  1  is  the  same  for  S  and  for  S’,  and 
therefore  processing  continues  identically  for  S  and  S’  after  that  point.  Thus, 
SELECT(S,c,p,i)  =  SELECT(S’,c,p,i). 

Case  2:  Bin  s  contains  more  than  one  resource  after  processing  S(1)...S(j-1)  on  c  and 
p,  or  else  c(choice(S,c,p,j))  is  not  bin  s.  (That  is,  the  bin  selected  by  the  choice  made  at 
step!  is  not  bin  s.) 

Then  the  effect  of  the  pair  of  steps  j  and  i  + 1  is  the  same  for  both  S  and  S’:  a  resource 
is  removed  from  bin  s  and  a  resource  iS  removed  from  bin  c(choice(S,c,p,j)).  (When 
processing  S,  the  choice  from  c  occurs  first,  while  when  processing  S’,  the  explicit 
removal  from  s  occurs  first,  but  the  net  effect  is  the  same.)  Subsequent  processing  of  the 
two  scripts  will  be  identical,  and  once  again,  SELECT(S,c,p,i)  =  SELECT(S’,c,p,i). 

Case  3:  Bin  s  contains  exactly  one  resource  after  processing  S(1)...S(j-1),  and 
c(choice(S,c,p,j))  =  s.  (That  is,  the  bin  selected  by  the  choice  made  at  step  j  is  s.) 

In  this  case,  the  processing  of  S  uses  choices  from  c  at  both  steps  j  and  j  +  1,  because 
the  choice  of  s  at  step  j  removes  the  last  resource  from  bin  s,  and  so  a  choice  must  also  be 
made  at  step  j  +  1.  The  proce.ssing  of  S’  does  not  need  a  choice  at  step  j,  although  it  is 
forced  to  choose  by  the  X  at  step  j  +  1.  Thus,  in  both  cases,  step  j  removes  the  last 
resource  from  bin  s,  while  step  j  +  1  makes  a  choice  using  c.  Then  choice(S,c,p,j  +  1)  * 
choice(S’,c,p,j  +  1):  that  is,  the  same  entry  in  c  is  used  at  step  j  +  1  in  both  cases.  The 
combined  effect  of  steps  j  and  j  +  1  on  the  bins  is  the  same  for  the  two  scripts.  Subsequent 
processing  is  again  identical,  so  SELECT(S,c,p,i)  =  SELECT(S’,c,p,i)  for  bin  i  ^  s,  and 
SELECT(S,c.p,s)  =  SELECT(S',c,p,s)  +  1  >  SELECT(S’,c,p,s).l 

We  can  now  state  the  main  result  of  this  section. 

Corollary  7:  For  any  scripts  S  and  S’  such  that  S  dominates  S’,  and  for  any  choice 
sequence  c,  placement  function  p  and  bin  i, 

SELECT(S,c,p,i)  >  SELECT(S’,c,p,i). 

Proof;  Let  T  be  an  augmentation  of  S’  by  a  suffix  of  X’s,  such  that  T  is  reachable  from 
S.  Then  Lemma  6  implies  that  SELECT(S,c,p,i)  >  SELECT(T,c,p,i).  But  the  latter  term  is 
obviously  at  least  as  great  as  SELECT(S’,c,p,i).l 


4.5.  Expansior}S 

In  this  subsection,  we  show  that  the  number  of  "BUYER"  messages  sent  is  a  monotone 
nondocroasing  function  of  the  intorartival  time  of  the  arrival  distribution.  Wo  do  this  by  comparing 
particular  pairs  of  e.xccutions. 

For  n  (i  N,  let  [n]  denote  {1 . n}.  If  a  C  R  and  r  =  (v,L)|(:  j|^j,  is  a  request  pattern,  then  ar,  the 

expansion  of  r  by  a,  is  the  request  pattern  (V|.at|)j^ji^|.  That  is,  ar  represents  the  request  pattern  in 


25 


which  the  successive  requests  occur  at  the  same  locations  as  in  r,  but  in  which  the  times  are 
expanded  by  the  constant  factor  a. 

We  will  compare  executions  for  request  pattern  r  and  request  pattern  ar,  using  the  same  choice 
sequence.  We  require  a  technical  restriction,  just  to  avoid  the  complications  of  having  to  consider 
multiple  events  occurring  at  the  same  node  at  the  same  time,  in  either  execution.  A  requesi  pattern,  r, 
is  said  to  be  a-isolated  provided  that  no  two  requests  occur  in  r  at  the  same  time,  and  provided  that 
the  following  holds.  If  t.|  and  are  two  times  at  which  requests  arrive  in  r,  where  t.|  t2,  ar  d  if  k  is  an 

integer,  then  the  following  are  true:  (a)  -  tj  ^  2k,  and  (b)  ■  tg  (2/a)k. 

The  next  lemma  is  crucial  to  our  analysis.  Its  truth  was  first  observed  empirically,  and  tfien  proved 
analytically.  It  says  that  the  number  of  "BUYER"  messages  sent  during  an  execution  cannot  increase 
if  the  request  pattern  is  expanded  by  a  constant  which  is  greater  than  or  equal  to  1 . 

Lemma  8:  If  a  >  1,  and  r  is  a  isolated,  then  bnum.|.  p(r,C)  <  bnum.p  p(ar,C). 

Proof;  Fix  T,  p  and  C.  Let  bsent(r,e,t)  denote  the  number  of  "BUYER"  message.3  sent 
along  edge  e,  in  the  execution  for  r  (using  T,  p  and  C),  up  to  and  including  time  '.  Let 
brec(r,v,t)  denote  the  number  of  "BUYER"  messages  received  by  vertex  v,  in  the  execution 
for  r,  up  to  and  including  time  t.  Let  arec(r,e,t)  denote  the  number  of  "ARRIVAL" 
messages  received  along  edge  e,  in  the  execution  for  r,  up  to  and  including  time  t.  We  will 
show  the  following: 

Claim: 

bsent(r,e,t  +  hGight.j(e))  <  bsent(ar,e,at  +  height.|.(e))  for  all  r,  e,  and  t. 

This  is  a  stronger  claim  than  required  for  the  lemma,  since  it  shows  an  inequality  not 
only  for  the  total  number  of  "BUYER"  messages,  but  for  the  number  along  each  edge,  up 
to  corresponding  times. 

Fact  1:arec(r.e,t  >  hc'ight.^(e))  =  arec(ar,e,at  +  height.^(e)). 

This  is  so  because  the  number  of  requests  arriving  in  request  sequence  r  by  time  t  is 
the  same  as  the  number  arriving  in  request  sequence  ar  by  time  at,  and  messages  just  flow 
up  the  tree  at  a  steady  r.ite. 

The  rest  of  the  proof  proceeds  by  induction  on  height^(e),  starting  with  height|.(e)  = 
height  p  and  -vorking  downv.ard  towards  the  leaves. 

Base:  hoight.^(e)  height.^ 

In  this  case,  e's  upper  endpoint  is  root^.  The  actions  of  rootj.  are  completely 
determined  by  the  "AHitIVAL"  messages  it  receives,  which  are  the  same  at  corresponding 
times  in  the  two  executions,  by  Fact  1.  The  Claim  follows. 


Inductive  5dop:  height  j(t')  <  height.j. 
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Let  V  be  the  upper  endpoint  of  edge  e,  and  the  lower  endpoint  of  edge  e'. 

Fact  2:  brec(r,v,t  +  height-p(v))  <  brec(ar,v,at  +  height^(v)). 

This  is  so  because  brec(r,v,t  +  height^(v))  =  bsent(r,e',t  +  height^(v)  •  1)  because  all 
messages  take  exactly  time  1,  =  .bsent(r,e’,t  -  2  +  height^(e')),  <  bsent(ar,e',a(t-2)  + 
height-^(e'))  by  inductive  hypothesis,  =  bsent(ar,e’.a{t-2)  +  1  +  height^(v)),  » 
brec(ar,v,a(t-2)  +  2  +  height^(v))  <  brec(ar,v,at  +  height^(v)),  since  a{t-2)  +  2  <  at. 

Now  let  us  consider  the  situation  from  v’s  viewpoint.  Node  v  is  comparing  two 
executions,  the  first  for  r  and  the  second  for  ar.  All  v  sees  is  its  incoming  "ARRIVAL"  and 
"BUYER"  messages,  and  v  uses  the  same  choice  sequence  in  both  cases.  At 
corresponding  times  in  the  two  executions  (i.e.  t  +  height.j.(v)  in  the  first  execution  vs.  at  + 
height.|.(v)  in  the  second  execution),  the  same  number  of  "ARRIVAL"  messages  have  been 
received  along  each  edge  (by  Fact  1),  and  an  inequality  holds  for  the  number  of  "BUYER" 
messages  which  have  been  received  (by  Fact  2).  We  will  show  the  needed  inequality  for 
the  number  of  "BUYER "  messages  sent  by  v  along  each  edge,  up  to  corresponding  times. 

Fix  any  time  t.  We  compare  the  first  execution  up  to  time  t  +  height.|.(v)  with  the 
second  execution  up  to  time  at  +  height.j.(v).  We  claim  that  this  situation  is  modelled  by 
the  combinatorial  problem  presented  in  the  preceding  subsection.  First,  we  represent  v’s 
inputs  in  each  of  the  two  executions  by  a  script,  i.e.  a  sequence  of  X's  and  "bins"  the  latter 
of  which  are  identified  with  children  of  v.  An  X  models  the  arrival  of  a  "BUYER",  while  s  € 
bins  models  the  receipt  of  a  "ARRIVAL"  message  from  s.  (The  fact  that  r  is  a-isolated 
means  that  no  two  of  v's  inputs  occur  at  the  same  time  in  the  same  execution,  so  a  unique 
sequence  can  be  obtained  in  each  case.)  Let  S  and  S’  be  the  scripts  for  the  first  and 
second  executions,  respectively  (up  to  the  indicated  times).  The  claims  in  the  preceding 
paragraph  imply  that  S’  dominates  S. 

V/e  claim  that  the  processing  described  for  the  combinatorial  problem  models  the 
procesring  carried  out  at  v  during  execution  of  the  algorithm.  In  particular,  a  SELECT  of  a 
bin  s,  if  It  occurs,  models  the  sending  of  a  "BUYER"  message  to  s  and  associated 
I  eduction  of  v's  estimate  of  tlie  number  of  resources  remaining  in  s’s  subtree.  With  the 
given  correspondence  between  the  combinatorial  problem  and  the  executions,  the 
conclusion  of  Corollary  7  translates  immediately  into  the  Claim. I 

Lemma  9;  If  a  >  1,  and  r  is  a  isolated,  then  bnumT-  (r)  <  bnum,.  (ar). 

—  T,p'  —  r,p'  ' 

Proof:  By  Lemma  8,  taking  expectations.! 

Define  bnum^^(Hj)  to  bo  the  expected  value  of  bnum.^p(ar),  where  r  is  chosen  according  to  (p.^ 
and  f.  The  next  theorem  states  that  the  expected  number  of  "BUYER"  messages  is  a  monotone 
nondecreasing  function  of  the  interarrival  tune  of  the  request  distribution. 

Theorem  10:  (a)  if  a  >  1,  then  bnum^^(f)  <  bnumj.  (a,f).  (b)  If  0  <  a  <  b,  then 
bnum,  (a.f)  <  bnum,  (b,f). 

i,p  —  i,p'  ' 


Pioof:  (a)  If  a  request  sequence  is  chosen  according  to  q>j  and  f,  then  with  probc.bility 
1,  it  will  be  a-isolated.  The  result  then  follows  from  Lemma  9,  by  taking  expectations  over 
r. 

(b‘  Let  g  be  the  probability  density  function  defined  by  g(at)  =  f(t).  Since  b/a  >  I,  we 
can  aoply  Part  (a)  to  b/a  and  g,  obtaining  bnum^  (g)  <  bnum^  (b/a,g).  But  bnum .  (g) 

=  bnum-j.  p(a,f)  and  bnum^  p(b/a,g)  =  bnum^p(b,f),  yielding  the  result. I 

4.6.  Summary  of  Monotonicity  Analysis 

In  this  s€  ction,  we  have  made  the  following  restrictions,  repeated  here  for  convenience. 

Restriciion  1: 

T  is  a  weighted  tree,  and  the  nondeterministic  choice  step  in  Part  (2)  of  the  algorithm  uses  a  call  to 
random.j.g. 

Restric  :ion  2: 

Delivery  time  for  messages  is  always  exactly  1. 

Restric  ion  3: 

T  has  all  leaves  at  the  same  distance  from  the  root,  and  r  and  p  are  nonzero  only  at  leaves. 

The  major  results  we  have  proved  in  this  section  are  that  the  expected  response  time  is  closely 
approximated  by  the  expected  number  of  "BUYER"  messages  (Theorem  5)  and  that  the  expected 
number  of "  3UYER"  messages  is  a  monotonic  function  of  the  interarrival  time  (Theorem  lOi.  We  can 
combine  these  two  results,  obtaining  the  following: 

Th  eorem  1 1 ;  cost,  (f)  <  4  lim  n^ibnumT  (a,f). 

Proof:  Consider  what  happens  to  the  value  of  bnum.^p(a,f)  as  a  increases.  This  value 
is  monotonically  nondecreasing,  by  Theorem  10.  Also,  it  is  bounded  above,  because  no 
execution  causes  more  "BUYER"  messages  to  be  sent  on  any  edge  than  the  number  of 
resources  in'iin.lly  placed  below  that  edge.  Therefore,  the  limit  exists.  The  result  follows 
immediately  from  Theorems  5  and  10.1 

That  is,  the  expected  cost  of  the  algorithm  for  any  probability  function,  f,  is  bounded  in  terms  of  the 
limiting  case  of  the  algorithrri,  as  the  inlerarrival  time  approaches  infinity.  But  note  that  as  the 
interarrival  t  me  apirroaches  infinity,  the  algorithm  gravitates  towards  a  purely  sequential  algorithm  • 
one  m  which  each  request  gets  satisfied  before  the  next  one  arrives.  This  kind  of  sequential  algorithm 
is  am  'na.ble  to  analysis  of  a  more  traditional  kind,  the  subject  of  the' next  section  of  this  paper. 


It  seems  quite  surprising  that  the  sequential  case  is  the  worst  case.  Our  initial  expectation  was  that 
cases  where  considerable  interference  between  requests  occur  would  be  the  worst  case.  The 
monotonicity  theorem  indicates  that  that  is  not  so.  Of  course,  we  have  made  a  few  assumptions  in 
this  section,  most  significantly  the  equal  lengths  of  branches.  It  is  quite  likely  that  the  sequential  case 
will  not  be  the  worst  case  for  an  algorithm  using  more  general  tree  topologies.  The  analysis  in  this 
more  general  case  so  far  seems  quite  intractable,  however. 

5.  Sequential  Analysis 

In  this  section,  we  analyze  the  sequential  case  of  the  algorithm.  In  the  next  section,  we  combine 
the  results  into  an  upper  bound  for  the  entire  algorithm.  Once  again,  we  allow  arbitrary  weighted 
trees  T,  and  allow  r  and  p  to  be  nonzero  anywhere. 


5.1 .  A  Simplified  Problem  and  Algorithm 

The  sequential  case  of  the  algorithm  offers  considerable  simplification  over  the  concurrent  cases. 
Tfiere  is  no  interference  at  all,  since  each  request  arrives  after  previous  requests  have  been  satisfied. 
This  means  that  all  the  estimates  of  remaining  resources  are  completely  accurate.  In  fact,  the  result  is 
equivalent  to  that  of  an  algorithm  in  which  ail  information  is  known  globally. 

The  behavior  of  the  algorithm  in  the  sequential  case  can  be  modelled  by  repeated  calls  to  the 
following  sequential  program,  FIND.  The  program  takes  a  weighted  tree,  a  nonnull  placement,  and  a 
vertex  (the  vertex  at  which  a  request  occurs)  as  input,  and  returns  a  vertex  (the  vertex  at  which  the 
resource  to  be  granted  is  located). 

FIND(T,p,v) 

Case 

p(v)  >  0  :  return  v 

p(desCj(v))  =  0  :  return  FIND{T,p,parentf(v)) 

p(v)  =  0  and  p(desCf(v))  >  0: 

[S  :=  {w  €  ch i  1  dren j^( v) :  p{desCf(w))  >  0} 
return  F IND( T , p , random^  j) ] 
endcase 


Thus,  a  request  is  satisfied,  as  before,  in  the  smallest  containing  subtree  which  contains  a 

resource:  where  there  is  a  choice,  the  probability  function  is  used. 

Lemma  12:  If  p  is  nonnull  and  v  C  vertices.|.,  then  FIND(T,p,v)  eventually  halts  and 
returns  a  vertex,  v,  with  p(v)  >  0. 


Proof:  Straightforward.! 


We  next  want  to  prove  a  lemma  which  will  be  useful  in  the  later  analysis.  The  content  of  the  lemma  is 
as  follows.  Let  random-^,  denote  the  probability  function  which  returns  s  €  leaves^  with  probability 
Assume  T  is  a  weighted  tree  and  p  is  a  placement  which  is  nonzero  only  at  leaves.  Consider 
the  following  two  experiments. 

(1)  Call  FIND{T,p,root.|.),  and 

(2)  Call  FIND(T,p, random^). 

We  claim  that  the  "results"  of  these  two  experiments  are  the  same.  That  is,  for  each  w  C 
vertices^,  the  probability  that  experiment  (1)  returns  w  is  exactly  the  same  as  the  probability  that 
experiment  (2)  returns  w.  It  will  follow  from  the  next  lemma  that  this  is  so. 

Some  notation  is  helpful  here.  The  result  of  FIND  on  a  particular  set  of  arguments  can  be 
expressed  as  a  probability  distribution  of  vertices.  Let  denote  the  probability  distribution  of 
results  for  FIND(T,p,v).  That  is,  FIND(T,p,v)  returns  w  with  probability  a.|.  ^  ^(w). 

Lemma  13:  Let  T  be  a  weighted  tree,  p  a  nonnull  placement  for  T.  Assume  that  p  is 
nonzero  only  at  leaves.  Then  the  following  are  true. 

(a)  If  V  €  internal^,  then  =  I^echHdren^tv)  II'PT(descT(w))/.p^(desc.r(v))]  ^ 

(b)  If  V  €  vertices,,  then  a,  ^  ^  ([9T(w)/'PT(desc,(v))]  a^  ^  J. 

P  roof;  In  the  proof,  we  assume  T  and  p  are  fixed  and  write  in  place  of  a,  ^  etc. 

(a)  We  consider  cases. 

Case  1:  p(desc,(v))  =  0 

Then  since  the  algorithm  immediately  calls  FIND  on  parent, (v),  we  see  that  = 

,  ,  ..  Similarly,  for  all  children,  w,  of  v,  we  have  «  =  a  ....  Since 

parent^(v)  ^  i  ^  parent-|.(v) 

^-  /.'€chiidren,(v)  ['PT(desc,(w))/<f),(desc,(v))]  =  1,  the  result  follows. 

Case  2;  p(desc,(v))  >0 

Since  we  are  assuming  that  v  C  internal,,  we  know  that  p{v)  =  0.  Let  S  = 

{w  €  children, (v):  p(desc,(w))  >  0}.  Then  S  0.  The  third  case  in  the  algorithm 
holds,  and  we  have  that 

“v  =  ^w€sffT(desc,(w))/<p,(desc,(S)))«J.  Now, 

^wCch,idrcn.(v)  [[T,(desc,(w))/<p,{desc,{v))]  «J 


=  [[<p^(desc^(w))/9^(desc^(v))]  a  J 

+  ^w€chiidren^(v).s  [[9T(desc,(w))/<p^(desc^(v))]  a  J. 

By  the  remark  above,  this  sum  is  equal  to 
[<P-p(desc^(S))/<]p^(desc^(v))] 

+  ^w€chiidren^(v).sft<PT(desc^(w))/<f.^(desc^(v))]  aj. 

If  w  €  children^(v)  •  S,  we  know  that  p(desc^(w))  =  0,  so  that  =  “parent  (w)  “  “v' 

So  the  sum  above  is  equal  to  ^ 

[<p^(desc^(S))/9-p(desc^(v))]  +  3  [9^(desc^{w))/(p^(desc^{v))], 

=  [(p^(desc^(S))/(p-^(desc^(v))]  +  [(9^(desc^(v))  ■  (p^(desc^(S)))/  9^(desc^{v))], 

=  «v 

(b)  We  proceed  by  induction  on  the  height  of  v. 

Base:  v  6  leaves^ 

Then  the  only  w  in  desc^(v)  is  v  itself,  so  the  sum  on  the  right  is  just  [9j(v)/(py(v)]  a^, 

=  as  needed. 

Inductive  step;  v  €  internal.j. 

Then  [[9T.(desc^(w))/9^(desc.p(v))]  o  J,  by  part  (a), 

=  ^w€ch,)dren^(v)  [[<PT(desc,(w))/9^(desc.,(v))l  [[9^(s)/(rT(desc,(w))]  aj], 

by  inductive  hypothesis, 

“  ^wCchildfen^(v)  ^sCdesc^(w)  “s^’ 

=  ^s€desc.^(v)  [(<PT(s)/<rr(desc.r(''))]  "J- 
as  needed. I 

Part  (b)  of  this  lemma,  with  v  =  root^,  proves  equivalence  of  the  two  experiments  described  prior  to 
the  lemma. 
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We  can  restrict  attention  to  "request  sequences”  in  place  of  request  patterns,  in  the  sequential  case 
of  the  algorithm.  .Assume  that  T  is  a  weighted  tree,  and  p  is  a  placement  for  T.  A  request  sequence,  r, 
is  a  sequence  of  elements  of  vertices^,  representing  a  sequence  of  request  arrival  locations. 
Similarly,  a  resource  sequence,  r,  is  a  sequence  of  elements  of  vertices.|-,  representing  a  sequence  of 
resource  locations.  In  either  case,  let  length(r)  denote  the  number  of  elements  in  a  secjuence.  A 
resource  sequence,  s.  is  compatible  with  a  placement,  p,  provided  that  |s  Vv)|  <  p(v)  for  each  v  € 
vertices^.  (That  is,  the  resource  sequence  grants  at  most  the  number  of  resources  placed  at  each 
vertex.)  If  r  is  a  request  sequence  and  p  is  a  placement  with  total(p)  >  length(r),  then  a  matching  of  r 
and  p  is  a  p.iir  m  =  (r.s),  where  s  is  a  resource  sequence  compatible  with  p  and  length(r)  =  lengthfs). 
A  matching  describes  the  successive  locations  of  resources  which  are  used  to  satisfy  a  sequence  of 
requests. 

Next,  we  define  a  probabilistic  program  which  takes  as  inputs  a  request  sequence,  r,  and  a 
placement,  o,  with  total(p)  >  length(r),  and  returns  a  matching  of  r  and  p.  The  procedure  simply  uses 
FIND  repeatedly. 

MATCH(T,p,r) 

for  i  =  1  to  length(r)  do 
[s(i)  :=  riND(T,p,r(i)) 
p(s{  i))  :=  p(s(i))  -  1] 

Theorem  14;  Lot  r  be  a  request  sequence,  p  a  placement  with  total(p)  >  length(r). 

Then  MATCH(T,p,r)  will  eventLially  halt  and  return  a  matching  of  r  and  p. 

Proof:  Straightforv/ard.l 

This  algorithm  is  designed  to  behave  ex.ictly  as  the  sequential  case  of  our  general  algorithm. 

5.2.  Cost  Measures 

Let  d!st  Ju.  v)  denote  the  tree  distance  between  u  and  v.  If  m  =  (r.s)  is  a  matciiing,  then 
-()s-r  (cij  >  dusl,(f(i)..sfi))  TIurs.  the  "sequential  cost"  is  just  the  sum  of  the  tree  distances 

I  .p  I  I 

between  succcs.sive  rf  quest;',  and  their  corresponding  resources. 

If  r  is  a  ;(:(iuo;,t  S'-ruierv  with  length(r)  <  total(p).  then  define  seqcostj ^(r)  to  be  the  expected 
v.ilue  of  PO'icost|  ^(m).  whim  in  is  constiucted  using  MA rCH(T,p,r).  Let  seqcosly^  denote  the 
expected  v;  loo  of  -enrost,  (i),  where  r  is  of  length  total(p),  with  its  successive  locations  chosen 
independently  acrorduni  to  ip^. 

In  the  reni.iindi  T  of  fOi  -,  (  tion,  we  analy/'e  oOqcostj. 
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5.3.  A  Useful  Recurrence 

In  this  subsection,  we  present  a  solution  to  a  system  of  recurrence  equations.  This  solution  will  be 
useful  in  later  subsections. 

Let  c  €  R  .  Define  G^:  N  x  R  ^  *  R  ^  by  the  equations: 

G^(0,t)  =  0,  and  G^(k,t)  =  max{G^(k-1,t^)  +  G^fk-l.tg):  t^  +  t^^t)  +  ck-/t.  for k  >  1. 

Lemma  1 6:  For  all  k  >  0,  the  following  are  true: 

(a)  The  function  mapping  t  to  G^{k.t)  is  concave  downward  and  monotone 
nondecreasing. 

(b)  If  k  >  1,  then  Gj,(k,t)  =  2G^(k-1,t/2)  +  ckVt. 

Proof:  We  proceed  by  induction  on  k.  The  base,  k  =  0,  is  trivial.  For  the  inductive 
step,  letk  >  1.  If  t^  +  t^  <  t,  then  G^(k-1,t^)  +  G^^fk-l.t^)  <  2G^(k-1,(t^  +  t2)/2),  since  the 
inductive  hypothesis  states  that  the  function  mapping  t  to  G^(k-1,t)  is  concave.  This  is  in 
turn  <  2G^(k-1  ,t/2).  by  monotonicity.  Therefore,  G^{k,t)  =  2G^(k-1,t/2)  +  ckVt,  showing 
(b).  Since  each  term  is  concave  and  monotone,  the  sum  is  also,  showing  (a). I 

Theorem  1 6:  G^(k,t)  <  (3  V  2  +  4)c  V 1 2\ 

Proof;  By  Lemma  15,  G^(k,t)  =  0  if  k  =  0,  and  =  2G^(k-1,t/2)  +  ckVt  if  k  >  1. 
Expanding  this  recurrence,  we  see  that  G„(k,t) 

=  c[Z,^g  ^.,2*(k-i)yt/2‘]forallk>0, 

Let  X  =  1/  /2,  n  =  2*'.  Then 

G^(k,t)  =  (cv/tyn/V2)[Zj^o . 

=  {c/tVn//2)[(1  +  kx’'^’ •  (k+1)xV(1-x)^]. 

=  (c/tv'n)/(/2(1  - 1/72)^11  +  kx''*'-(k+1)x''], 

=  (cv't7n)(3/2  +  4)[1  +  kx'^*’  •(k+1)x''l, 

=  (c  y  t  y  n)(3y  2  +  4)(1  +  (kx  •  k  •  1)x''], 

<  (c  y  t  y  n)(3  y  2  +  4),  since  kx  •  k  -  1  <  0, 

=  (cyt2'')(3y2  +  4).i 
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5.4.  Recursive  Analysis 

Now,  we  require  Restriction  3  and  a  new  assumption,  Restriction  4.  These  are  to  remain  in  force 
for  the  remainder  of  Section  5. 

Restriction  3: 

T  has  all  leaves  at  the  same  distance  from  the  root,  and  r  and  p  are  nonzero  only  at  leaves. 
Restriction  4: 


T  is  a  complete  binary  tree. 

If  T  is  a  jveighted  tree,  then  a  weighted  subtree,  T,  of  T  consists  of  a  subtree  of  T  together  with  a 
probability  function.  <p^.,  given  by  <Pj,(v)  =  (p.j.(v)/<j().|.(vertices.j-,).  That  is,  the  weights  of  the  subtrees 
are  just  nor  nalized  restrictions  of  the  weights  of  T.  If  T  is  a  weighted  binary  tree,  let  leltj  and  right j 
denote  the  designated  weighted  subtrees  of  T. 


If  T  has  height  at  least  1,  then  let  T^  and  denote  left.p  and  right.j.,  respectively.  Let  p^  and  Pg 

denote  p|T^  and  pjT^,  respectively.  If  r  is  a  request  sequence,  let  overflow ^ ^(r)  denote 

||r  ’(vortices,.  )|  •  p(vertices,  )|,  the  difference  between  the  number  of  requests  that  arrive  in  the  left 
11 

subtree  and  the  number  of  resources  placed  there.  Let  overflow^ ^  denote  the  expected  value  of 
overflow,,  ^fr),  where  r  is  a  sequence  of  length  total(p)  chosen  using  rpj. 


The  follov;ing  is  a  key  lemma  which  provides  a  recursive  breakdown  for  the  sequenti.al  cost.  It 
says  that  the  expected  cost  of  matching  breaks  down  into  costs  of  matching  vathin  the  twc  subtrees, 


plus  a  charge  for  the  requests  that  overflow  into  the  opposite  subtree. 

Lemma  1  7:  Assume  height,.  >  1.  Then  seqcost.,p  <  seqeost.,. 
2  tieight,  overflow.,  p. 


Pi 


+  seqeost.,  + 
'2'P2 


Proof:  For  any  particular  request  sequence,  r,  there  is  some  particular  number, 
overflovy,  (r).  of  requests  that  do  not  get  satisfied  within  their  own  subtree,  but  rather 
overflow  into  the  opposite  subtree  to  find  a  resource.  To  be  specific,  assume  that  it 's  the 
left  subtree  from  which  any  excess  requests  overflow.  Lot  r,  be  the  subsequence  of  r 
cons'  . ting  of  requests  arriving  in  T,,  truncated  to  length  p(T,).  Let  r^  be  the  subsequence 
of  r  consisting  of  r«  quests  arriving  in  T^.  Recall  that  seqeost, p  is  the  expectation  of  the 
search  cost  for  r  nongh  requests  to  exhaust  all  resources  present. 

Before  the  time  at  which  the  left  subtree  overflows,  the  algorithm  MATCH{T,p,r)  runs 
exactly  like  MATCH(T,.p,,r,)  within  the  left  subtree.  Requests  originating  wittrn  T, 
become  matched  to  exactly  the  same  resources  in  both  executions. 


Wc  now  con  irlor  the  right  subtree.  Requests  which  originate  within  the  right  subtree 
are  handled  in  th.'  algorithm  MATCH(T.p,r)  exactly  as  they  are  in  the  algorithm  for  T,,  and 
p  ,  Kovac.'er,  there  aie  also  overflow  requests  from  the  left  subtree,  which  enter  T^  at  its 
roof  rather  th.in  .at  its  leaves.  By  Lemma  13.  whenever  such  a  request  arrives,  its 
probability  of  being  matched  to  my  particular  resource  is  exactly  the  same  as  if  the 
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request  had  entered  at  a  random  leaf  of  T^.  All  the  requests  remain  independent,  and 
these  additional  requests  are  just  enough  to  exhaust  the  resources  in  T^. 

Now  assume  that  the  request  sequences,  r,  are  chosen  at  random  according  to 
They  result  in  subsequences  r^  and  r^  which  are  chosen  at  random  according  to  and 

9,  ,  respectively.  ’ 

2 

We  claim  that  seqcost.,.  p,  the  expected  cost  taken  over  all  r,  is  bounded  by 

seqcost.|.  p  ,  the  expected  cost  taken  over  random  sequences  in  the  left  subtree, 

+  seqcost,  ,  the  expected  cost  taken  over  random  sequences  in  the  right  subtree, 

ij.Pg 

+  2  height.|.  overflow.|.  p. 

The  third  term  allows  for  the  expected  overflow  of  requests,  and  assigns  them  the 
maximum  cost,  2  height.^. 

Consider  the  first  term,  (The  second  term  is  analogous.)  The  first  term  allows  for  the 
expected  cost  incurred  by  an  execution  of  MATCH  on  a  random  sequence  of  requests 
within  T^. 

In  case  T^  has  its  resources  exhausted  by  requests  originating  within  that  subtree,  this 
term  measures  exactly  the  expected  cost  for  the  matching  of  the  first  requests  in  T,|  to  all 
the  resources  in  T^.  This  term  ignores  the  cost  incurred  by  any  excess  requests 
originating  within  T^  which  do  not  get  matched  within  T,|.  However,  that  is  not  a  problem, 
since  those  costs  are  counted  by  the  third  term. 

In  case  T^  does  not  have  its  resources  exhausted  by  requests  originating  within  T^,  this 
term  13  actually  greater  than  needed  to  measure  the  expected  cost  of  matching  all  requests 
originating  in  Tj',  in  fact,  it  is  enough  to  measure  this  cost  of  matching  these  requests, 
intersper.sod  with  enough  other  random  requests  (arriving  at  the  leaves  of  T^)  to  use  up  all 
the  resources  in  T^.  W.i  have  already  argued  that  these  requests  behave  as  if  they  were 
interspersed  with  other  random  requests,  because  requests  arriving  at  the  root  match  in 
the  same  way  as  if  they  arrived  at  random  leaves.  In  this  case,  the  first  term  does  not 
account  for  the  cost  of  matching  those  requests  which  enter  T^  at  its  root.  However,  that  is 
not  a  problem  since  that  cost  is  covered  by  the  third  term. I 


5.5.  d-Fairness 

We  need  to  make  another  restriction  on  the  algorithm,  for  the  purpose  of  analysis.  In  particular,  it 
is  reasonable  that  the  behavior  of  the  algorithm  should  be  best  when  the  resources  are  distributed  in 
the  tree  in  seme  relationship  with  the  probability  distribution  governing  arrival  of  requests.  (The  paper 
[Fischer,  Griffeth,  Guibas  and  Lynch  (1981)]  considers  optimal  placements  of  resources  in  a  tree 
network.) 


For  d  €  R*.  we  say  that  a  placement,  p,  is  d-tair  provided  that  the  following  is  true.  Let  u,v  € 
verticeSp  where  u  €  desc.|-(v).  Let  =  <pp(desCp(v))  and  <p^  =  <pp(desCp(u)),  Then  |p(desc^(u))  - 
(<P2^*^'i)P(^®sc.|.(v))|  <  d  y  (<)P2/<p,)p(desCp(v)).  That  is.  for  each  subtree,  the  number  of  resources 
placed  in  each  of  its  subtrees  is  approximately  proportional  to  the  probability  of  arrivals  in  that 
subtree,  and  the  difference  is  proportional  to  the  "standard  deviation".  For  any  T,  if  t  and  d  are 
sufficiently  large,  then  techniques  in  [Fischer,  Griffeth,  Guibas  and  Lynch  (1981)]  can  be  used  to 
produce  d  fair  placements  of  t  resources  in  T. 

From  now  on  in  Section  5,  we  assume  the  following. 

Restriction  5: 

p  is  d  fair  (for  some  arbitrary  but  fixed  d). 

The  next  lemma  says  that  restrictions  of  d-fair  placements  are  also  d  fair. 

Lemma  18;  Let  p  be  a  d  fair  placement  for  T.  Let  T’  be  any  subtree  of  T,  and  p’  = 
p|vertices.|^.,  the  restriction  of  p  to  vertices^,.  Then  p’  is  d-fair  for  T’. 

Proof:  Let  u,  v  £  vertices^,  with  u  €  desCp(v).  Let  9.|  =  (jp.p(desCp(v)),  <p^  = 
()f).^(desc^(u)),  <jp',  =  (pp(desc^,(v)),  and  <p'^  =  (pp(desCp,(u)).  Then  = 

(p,/(f'^(vertices^,)  and  <p'^  =  (p2/9T^(verticeSp),  by  definition.  Therefore,  (p'(2)/(p’(1)  = 

V'Pr 

Note  that  p’(desc^.(u))  =  p(desc^(u)),  and  p’(desCp,(v))  =  p(desc.^(v)).  Thus, 
|p'(dosCp(u))  •  p’(desCp(v))| 

-  |p(desc^(u))  •  (’P2''(p,)  P(desc^(v))|, 

<d  y  (()32/')P,)P(desc^(v))  since  p  is  d-fair, 

=  d  y  (<j)'2/<p',)p’(desc.^,(v)),  as  needed.l 


The  final  lemma  of  this  subsection  bounds  the  expected  overflow  for  d-fair  placements. 
Lemma  19:  overflow.|.p  <  (6  -r-  d)  y  <p.|.(verticeSp  )total(p). 

Proof:  ||r  '(vertices^  )|  -  p(vertices.p  )[  <  ||r ’(vertices^  )|  -  (p^(vertices.|.  )  total(p)|  + 

|(f'j(vertices,  )  total(p)l  •  |)(vertices^  )|.  ’  ^  ^ 

'  1  1 

The  expected  value  of  the  first  of  these  quantities,  taken  over  r,  is  bounded  by 
6  y<)f'.j.( vertices^  )  total(p),  using  Lemma  3.1.5  of  [Fischer,  Griffeth,  Guibas  and  Lynch 
(1901)].  ' 

The  second  quantity  is  bounded  by  d  y  <pp(vertices^  )  total(p),  since  p  is  d-fair.l 
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5.6.  Sequential  Analysis 

Let  sizBj.  denote  the  number  of  vertices  in  T.  We  now  give  the  main  result  of  Section  5,  that 
seqcost.^  p  is  0(  ^  size.^  total(p)).  This  says,  for  example,  that  if  total(p)  is  proportional  to  size.^,  then 
the  total  cost  is  proportional  to  total(p).  This  implies  that  the  average  cost  per  request  is  just  a 
constant,  independent  of  the  size  of  the  network. 

In  order  to  prove  this  theorem,  we  require  the  following  restrictions,  repeated  here  for 
convenience. 

Restriction  3: 

T  has  all  leaves  at  the  same  distance  from  the  root,  and  r  and  p  are  nonzero  only  at  leaves. 

Restriction  4: 

T  is  a  complete  binary  tree. 

Theorem  20;  seqcost.|.p  is  0(7  size^  total(p)). 

(More  specifically,  scqcost.|.p  <  (3 V  2  +  4){{6  +  d)-/  2)V  total(p).) 

Proof:  By  Theorem  16,  it  suffices  to  show  that  seqcost.^p  <  Gp(height.p,total(p)), 
where  c  -  (6  +  d)  V  2. 

We  proceed  by  induction  on  height.|.. 

Base:  height.^  =  0 

Then  T  has  a  single  vertex,  and  seqcost.|.  p  =  0.  The  inequality  is  immediate. 

Inductive  Step:  height.j.  >  1 

Then  seqcostr  <  seqcost,  +  seqcostT  +  2  height,  overflow,  ,  by  Lemma 
r.p  —  ^  ^  Tj.Pg  T  T,p  ' 

'  '  ! 

<  seqcost^  +  seqcostj.  +  2  height.^  (6  +  d)V  <p.j.(vertices.^  )total(p),  by  Lemma 

19,  r  1  2'  2  1 

A  similar  inequality  holds  for  T^  in  place  of  T^  within  the  square  root.  Since  at  least  one 
of  (j'^(vertice3.p  ),  (p^(vertices^  )  is  no  more  than  1/2,  it  follows  that  seqcost^  p  < 

seqcost,  +  seqcost,  +  2heightr  (6  +  d) /(1/2Hotal(p), 

',.P,  >2^2  ' 

=  seqcost,  +  seqcost,  +  (6  +  d)/2height,  ytotal(p). 

',.p,  '2^2 
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By  Lemma  18,  we  can  apply  the  inductive  hypothesis,  which  implies  that  the  right  hand 
side  of  this  inequality  is  at  most  equal  to 

G^(height^-  l,total(p^))  +  G^(height-j.  ■  l.totaKp^))  +  (6  +  d) ■/ 2  height,^  7total(p). 

The  definition  of  G^  implies  that  this  latter  expression  is  at  most  equal  to 
G^(height-|.,total(p)),  as  needed, I 


6.  The  Final  Analysis 

In  this  section,  we  combine  the  monotonicity  analysis  and  the  sequential  analysis,  to  obtain  an 
upper  bound  for  the  expected  cost  for  the  algorithm. 


6. 1 .  Relationship  Between  the  Costs 

Now  we  require  the  following  restrictions. 

Rest  riction  1 ; 

T  IS  a  weighted  tree,  and  the  nondeterministic  choice  step  in  Part  (2)  of  the  algorithm  uses  a  call  to 
random^  g. 

Restriction  2: 

Delivery  ti.mo  for  messages  is  always  exactly  1. 

With  these  restrictions,  there  is  a  close  relationship  between  the  costs  of  our  general  algorithm 
and  the  cost  of  the  sequential  algorithm  MATCH, 

Lemma  21:5oqcost^p  =  lim^  QQ(bnum.j-p(a,f)  +  vbnum.|.  p(a,f)), 

P,-oof:  There  is  an  absolute  upper  bound  on  the  time  for  our  algorithm  to  satisfy  a 
single  recuiest.  in  the  absence  of  concurrent  requests.  Thus,  as  a  increases,  the 
prob.nhility  that  there  are  any  concurrent  requests  approaches  0. 

Th<'r>Morf?.  the  limiting  case  of  the  general  algorithm  behaves  like  MATCH.  There  is  no 
backtracking,  so  the  total  search  time  just  reduces  to  the  sum  of  the  distances  from 
requests  to  the  resources  which  satisfy  them.  This  sum  is  just  the  total  number  of  buyer 
anri  vTtual  buyer  messages. I 


Now  let  us  add  one  more  rootriction; 

Restriction  3: 

T  has  all  Icavi.s  at  the  same  distance  from  the  root,  and  r  and  p  are  non^ero  only  at  leaves. 


With  this  added  restriction,  we  can  prove  a  variant  of  the  preceding  lemma. 
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Lemma  22:seqcost^p  =  2lim^_^  o^Cbnum^ p(a,f)). 

Proof:  By  Lemm-c  ?1  and  4.1 

This  immediately  implies  the  following  bound  on  the  cost  of  the  general  algorithm. 
Lemma  23:  cost.^  ^(f)  <  2  seqcost.j. 

Proof:  By  Theorem  1 1  and  Lemma  22.1 


6.2.  The  Main  Theorem 

Now  we  are  ready  to  present  the  main  result,  the  upper  bound  for  the  expected  cost  for  the 
general  algorithm  In  order  to  apply  the  results  of  both  the  monotonicity  analysis  and  the  sequential 
analysis,  we  must  assume  the  restrictions  made  for  both  cases.  More  specifically,  we  assume  all  of 
Restrictions  1-5: 

Restriction  1: 

T  is  a  weighted  tree,  and  the  nondeterministic  choice  step  in  Part  (2)  of  the  algorithm  uses  a  call  to 
random.pg. 

Restriction  2: 

Delivery  time  for  messages  is  always  exactly  1. 

Restriction  3: 

T  has  all  leaves  at  the  same  distance  from  the  root,  and  r  and  p  are  nonzero  only  at  leaves. 

Restriction  4: 

T  is  a  complete  binary  tree. 

Restriction  5: 

p  is  d  fair  (for  some  arbitrary  but  fixed  d). 

Theorem  24:  Let  f  be  a  probability  function.  Then  cost.|.  ^(f)  is  0(7  size.j.  total(p)). 

(More  specifically,  cost^  p(f)  <  2(37  2  +  4)((6  +  d)7  2)7  total(p).) 

Proof:  By  Lemma  23  and  Theorem  20.1 

In  particular,  provided  that  total(p)  is  proportional  to  size.^,  the  expected  average  time  taken  by 
this  alyorilhm  to  .satisfy  a  single  request  is  constant,  independent  of  the  size  of  the  network. 
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Remark:  It  is  possible  to  prove  a  variant  of  Theorem  24,  for  the  case  in  which  the  placemt  i  p  is 
chosen  at  random  using  (just  as  the  request  locations  are  chosen),  rather  than  being  d-fair.  We 
sketch  the  ideas  briefly. 

First,  we  must  extend  the  cost  definitions  to  include  expectations  taken  over  placements  of  a 
particular  length.  Thus,  we  define  cost.^  ,(f)  to  be  the  expected  value  of  cost.^  p(f)  for  p  with  totalfp)  = 
t.  Analogous  definitions  are  made  for  seqcost.^j  and  overflow.j.  j.  Lemma  23  then  implies  that 
cost.^  ,(f)  <  2  seqcost.|.  It  is  also  easy  to  see  that  overflow.|. ,  <  d  7  <p.|.(vertices.^  )  t,  for  some 
constant  d. 


Next,  we  prove  a  consequence  of  Lemma  17  which  says  that  seqcost^  ^  <  Exp^^  ^  ^  ^  t  -  t 

(seqcosU  +  sepcost^.  ,  )  +  2  height,  overflow,,.  (Here,  the  expectation  is  taken  over  pairs 
r  1  '2' 2  '•* 

which  are  obtained  by  using  <p.j.  to  assign  resources  to  T.|  or  Tg.)  This  obviously  implies  that 
seqcost^  j  <  max^,  .t^iwitht  +  ,t  seqcost.^  ,  )  +  2  height^  overflow.p  j. 

Now  we  prove  a  variant  of  Theorem  20  which  says  that  seqcost.pj  is  0(7  size^  t).  More 
specifically,  we  show  that  seqcost.j.,  <  GJheight.p,t),  where  c  =  d  7  2.  This  is  easily  done  by 
induction  as  before,  using  the  new  lemmas  just  described. 

Combining  these  results,  we  see  that  cost.j.  ^(f)  is  0(  7  size.^  t). 


7.  Remaining  Questions 

There  are  several  directions  for  remaining  research. 

First,  we  v;ould  like  to  extend  the  analysis  of  the  general  algorithm.  We  would  like  to  loosen  our 
restrictions  on  tree  shape,  message  delivery  time,  and  locations  for  resources  and  requests.  If  we  do 
this,  IS  it  possible  to  carry  out  an  analysis  similar  to  the  one  in  this  paper?  In  particular,  can  the 

i  .irrent  ca.se-;  of  tfie  algorithm  still  be  bounded  in  some  way  in  terms  of  the  sequential  case? 

'vVe  would  also  like  to  extend  the  analysis  of  the  sequential  algorithm,  MATCH.  Here,  we  would 
'  ■  '  >0  'oosKfi  restrictions  on  tree  sh.ape  and  on  locations  for  resources  and  requests. 

I  here  are  some  apparent  improvements  in  the  algorithm,  for  example  adjusting  the  probabilities 
f  'r  !hi>  choice  among  children  in  response  to  knowledge  of  the  number  of  re.sources  remaining  in 
■  ■  u  n  subtree,  vv:  ile  this  seems  like  an  improvement,  the  resulting  algorithm  seems  harder  to  analyze 
( ,111^0  (he  recursive  decomposition  doesn't  appear  to  work).  Can  any  simple  modifications  be  shown 
to  be  improvements? 

Wt'  would  like  to  compare  the  performance  of  this  algorithm  to  that  of  alternative  algorithms  which 
;,nlv(>  the  same  problem.  We  have  alretidy  observe<l  that  this  algorithm  performs  much  better  than  the 
Cl  lit/  ili.-ed  .alfjorithm.  which  locates  all  resources  at  the  root.  How  does  it  compare  to  algorithms 


40 


which  allow  requests  to  search  for  resources  in  parallel  rather  than  sequentially?  What  about 
algorithms  which  rebalance  resources?  Are  there  other  interesting  ideas  for  algorithms  for  this 
problem? 

Finally,  the  general  analysis  strategy  is  quite  attractive.  Proving  a  monotonicity  result  which 
bounds  the  concurrent  cases  of  an  algorithm  in  terms  of  the  sequential  case,  and  then  analyzing  the 
sequential  case  by  traditional  techniques,  appears  quite  tractable.  The  use  of  this  strategy  for  our 
algorithm  appears  to  depend  on  many  special  properties  of  the  algorithm  and  on  restrictions  on  the 
execution.  Is  the  strategy  more  generally  useful?  For  what  type  of  algorithms  can  it  be  used? 

8.  References 

Fischer,  M.,  Griffeth,  N.,  Guibas,  L.,  and  Lynch,  N.  (1981),  Optimal  placement  of  identical 
resources  in  a  distributed  network,  in  "Proc.  2nd  International  Conference  on  Distributed 
Computing”.  Also,  to  appear  in  Inf  or.  and  Control. 

Guibas,  L.J.,  and  Liang,  F.  M.  (1982),  Systolic  stacks,  queues,  and  counters", 
in  "Proc.  2nd  MIT  VLSI  Conference". 


OFFICIAL  DIETRIBUTTO:.  LIST 

1985 

Director 

Information  Processing  Techniques  Office 

Defense  Advanced  Research  Projects  Agency 

1400  Wilson  Boulevard 

Arlington,  VA  22209 

2 

Copies 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  VA  22217 

Attn:  Dr.  R.  Grafton,  Code  433 

2 

Copies 

Director,  Code  2627 

Naval  Research  Laboratory 

Vv’ashing ton ,  DC  2  0  375 

6 

Copies 

Defense  Technical  Information  Center 

Cameron  Station 

Alexandria,  VA  22314 

12 

Copies 

National  Science  Foundation 

Office  cf  Conputinc  Activities 

1800  G.  Street,  N.K. 

Nash  me  ton,  DC  205  5  0 

Attn:  Program.  Director 

2 

Copies 

Dr.  L.B.  Royce ,  Code  38 

Head,  Research  Department 

Naval  Weapons  Center 

China  Lake,  CA  *'^  35  5  5 

1 

Copy 

Dr.  G.  Hopper,  USNR 

NAVDAC-OOH 

Department  of  t.he  Nax’y 

Washington,  DC  20374 

1 

Copy 

